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Dihadron azimuthal correlations containing a high transverse momentum (pr) trigger particle are 
sensitive to the properties of the nuclear medium created at RHIC through the strong interactions 
occurring between the traversing parton and the medium, i.e. jet-quenching. Previous measurements 
revealed a strong modification to dihadron azimuthal correlations in Au+Au collisions with respect 
to p+p and d+Au collisions. The modification increases with the collision centrality, suggesting 
a path- length dependence to the jet-quenching effect. This paper reports STAR measurements of 
dihadron azimuthal correlations in medium central (20-60%) Au+Au collisions at ^/s NN = 200 GeV 
as a function of the trigger particle's azimuthal angle relative to the event plane, <f> a = \4>t — V'epI- 
The azimuthal correlation is studied as a function of both trigger and associated particle pr- The 
subtractions of the combinatorial background and anisotropic flow, assuming Zero Yield At Mini- 
mum (zyam), are described. The away-side correlation is strongly modified, and the modification 
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varies with <f> s , which is expected to be related to the path-length that the away-side parton tra- 
verses. The pseudo-rapidity (At]) dependence of the near-side correlation, sensitive to long range 
At] correlations (the ridge), is also investigated. The ridge and jet-like components of the near-side 
correlation are studied as a function of 4> s . The ridge appears to drop with increasing <p s while 
the jet-like component remains approximately constant. These results, with extensive systematic 
studies of the dihadron correlations as a function of <j> s , trigger and associated particle pr, and 
the pseudo-rapidity range Ar], should provide stringent inputs to help understand the underlying 
physics mechanisms of jet-medium interactions in high energy nuclear collisions. 

PACS numbers: 25.75.-q, 25.75.Dw 



I. INTRODUCTION 

Collisions at the Relativistic Heavy Ion Collider 
(RHIC) of Brookhaven National Laboratory have created 
a medium with properties that resemble a nearly perfect 
liquid of strongly interacting quarks and gluons [H-Q- 
This conclusion is based upon two pillars of evidence: (i) 
the strong elliptic flow and (ii) jet-quenching- suppres- 
sion of high transverse momentum (px) single hadron 
yield and dihadron correlation in heavy-ion collisions rel- 
ative to elementary p+p interactions. While suppres- 
sion of high pt single hadron yields has limited sen- 
sitivity to the medium core, dihadron correlation mea- 
surements provide richer and more valuable information 
about the properties of the created medium [||. There 
are several key observations that can be made from di- 
hadron correlations with a high px trigger particle, (i) 
The correlated hadron yield at high px, while not much 
changed on the near side of the trigger particle (where 
azimuth difference between correlated and trigger par- 
ticles \A(f>\ < ir/2), is strongly suppressed on the away 
side (where \A(j)\ > ir/2) Q. This lends strong support 
to the partonic energy loss picture 043- The cor- 
related hadron yields at low px are strongly enhanced 
on both the near and away side [To| . In particular, 
the near-side enhancement is tied to lon g-ra nge corre- 
lations in pseudo-rapidity - the ridge [ToMT^ I . (hi) The 
away-side correlation broadens from peripheral to cen- 
tral collisions, and exhibits double peaks for select trig- 
ger and associated particle px ranges [1(1 [H, EH ■ The 
double peak structure opens up a new opportunity to 
study the underlying physics mechanisms for partonic en- 
ergy loss, such as gluon radiation [l5|, [l6j or Mach-cone 
shock- wave excitation (T7rl20| . Recent three- particle jet- 
like correlation studies indicate that the conical emission 
pattern of correlated hadrons is characteristic of Mach- 
cone shock- waves [2l| . (iv) The away-side associated par- 
ticles are partially equilibrated with the bulk medium 
in mid-central to central collisions, and a higher degree 
of equilibration is observed for particles which are more 
aligned back-to-back with the trigger particles [Io|, EH ■ 
This observation underscores the connection between the 
medium's path-length and partonic energy loss. 

We study the path-length dependence of partonic en- 
ergy loss in detail in non-central collisions where the over- 
lap region between the two colliding nuclei is anisotropic: 
the size in the reaction-plane direction is shorter than 



that perpendicular to it. The reaction plane (RP) is de- 
fined by the beam direction and the line connecting the 
centers of two colliding nuclei. It can be estimated in 
non-central collisions by determining the azimuthal an- 
gle with the highest particle emission probability, using 
the fact that the particles have an elliptic emission pat- 
tern. The estimated angle is called the event plane (EP) 
to emphasize that it is an experimental estimate of the 
reaction plane with finite resolution. By selecting the 
trigger particle direction with respect to the event plane, 
4> s = \4>t — V'epI (where 4>t is the trigger particle azimuth 
and tpEP is the event-plane azimuth), we effectively se- 
lect different path-lengths through the medium that the 
away-side parton traverses, providing differential infor- 
mation unavailable to inclusive jet-correlation measure- 
ments. 

Previously, the Solenoidal Tracker at RHIC (STAR) 
experiment has performed an exploratory measurement 
of azimuthal correlations at high px with trigger particles 
in-plane (4> s < 7r/4) and out-of-plane (4> s > 7r/4) using 
non-central 20-60% Au+Au collisions [22|]. The results 
hinted that the away-side correlation with out-of-plane 
trigger particles is more strongly suppressed than that 
with in-plane trigger particles. In this paper, we extend 
those measurements to finer bins in cj) s and to lower as- 
sociated and trigger px ranges [23[ . We also present in- 
clusive jet-correlation results from minimum bias d+Au 
collisions as a reference to the Au+Au data. 

We further study the ridge as a function of 4> s , and 
investigate the systematics of the ridge in an attempt to 
identify the underlying physics mechanism for the forma- 
tion of the ridge. 

The paper is organized as follows. In Sec. |TT] we de- 
scribe in detail our data analysis of dihadron correlations 
relative to the event plane. In Sec. IIIII we discuss our 
extensive studies of the systematic uncertainties of our 
results. In Sec. lIVI we report and discuss our results of di- 
hadron correlations relative to the event plane. We finally 
conclude in Sec. [V] We present all raw and background- 
subtracted dihadron correlation functions relative to the 
event plane in Appendix IA1 



II. EXPERIMENT AND DATA ANALYSIS 

The data used in this analysis were taken by the STAR 
experiment [24| at RHIC at the nucleon-nucleon center 
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of mass energy of ^/s NN = 200 GeV. The minimum- 
bias Au+Au data were from Run IV in 2004 at RHIC. 
The reference minimum-bias d+Au data used for com- 
parison were from Run III in 2003. The minimum-bias 
triggers for Au+Au and d+Au collisions were provided 
by the Central Trigger Barrel [25| and the Zero Degree 
Calorimeters [2a |. 

The details of the STAR experiment can be found in 
Ref . [24| • The main detector used for this analysis is the 
Time Projection Chamber (TPC) (HSI]. The TPC is 
surrounded by a solenoidal magnet providing a nearly 
uniform magnetic field of 0.5 Tesla along the beam direc- 
tion. Particle tracks are reconstructed in the TPC. The 
primary event vertex was fit using reconstructed tracks 
which pass certain quality cuts. Events with a primary 
vertex within ±30 cm of the geometric center of the TPC 
along the beam axis are used in the analysis. With this 
range of primary vertex position, the TPC has good ac- 
ceptance within the pseudo-rapidity region of |ry| < 1.1. 
The Au+Au collision centrality is defined according to 
the measured charged hadron multiplicity in the TPC 
within < 0.5 (reference multiplicity) [2{|. We choose 
the 20-60% centrality data for our analysis where good 
event-plane resolution is achieved (see later). 

In our analysis, only tracks that extrapolate to within 
2 cm of the primary vertex are used. Tracks are required 
to be reconstructed with at least 20 out of a maximum of 
45 hits in the TPC. The ratio of the number of hits used 
in track reconstruction to the number of possible hits 
is required to be greater than 0.51 to eliminate multiple 
track segments being reconstructed from a single particle 
trajectory. The same event and track cuts are applied to 
particle tracks used for event-plane reconstruction and 
for the subsequent correlation analysis. Particle tracks 
within 1 77 1 < 1 are used in the correlation analysis. 

High pt particles are selected as triggers off-line to 
perform the correlation analyses. We select high pt trig- 
ger particles within the pj) ranges of 3-4 GeV/c and 
4-6 GeV/c. A total of 4.4 million Au+Au events with 
centrality ranging from 20-60% are used in this analy- 
sis. From the event sample we find 2.1 million trigger 
particles with pr values ranging between 3-4 GeV/c and 
0.36 million trigger particles with pt values between 4- 
6 GeV/c. Associated particles, i.e. all particles in the 
event including those correlated with the trigger parti- 
cles, are grouped into the p^ ' ranges of 0.15-0.5, 0.5-1.0, 

1.0-1.5, 1.5-2.0, 2.0-3.0 (or 2.0-4.0) GeV/c. The low pfi 
cut-off of 0.15 GeV/c was imposed by the magnetic field 
strength and the TPC aperture. The azimuthal correla- 
tion functions in A</> (azimuthal angle difference between 
associated particle and trigger particle) are analyzed sep- 
arately for trigger particles at different azimuthal angles 
((j) s ) relative to the event plane. 

The associated particle yields are corrected for single- 
particle track reconstruction efficiency which is obtained 
from embedding simulated tracks into real events [30j . 
It depends on both centrality and pr- The efficiency is 



found to be insensitive to 77 and is therefore averaged over 
77. The </>-dependent part of the acceptance and track re- 
construction efficiency are corrected for both the trigger 
and associated particle yields. This (^-dependent correc- 
tion is obtained from the inverse of the single-particle 4> 
distribution whose average is normalized to unity. Cor- 
rection for the ^-dependent efficiencies for both trigger 
and associated particles removes the majority of the non- 
uniformity caused by the TPC sector boundaries. The 
remaining non-uniformity in A<j) is corrected by using an 
event-mixing technique, where the trigger particle from 
one event is paired with associated particles from an- 
other event within the same centrality bin [10j . The two- 
particle acceptance in A77 (pseudo-rapidity difference be- 
tween associated particle and trigger particle), which is 
approximately triangle-shaped in A77, is not corrected to 
be consistent with earlier publications The correla- 
tion function is normalized by the corrected number of 
trigger particles in its corresponding (f> s bin. The cen- 
trality and pr dependent aspects of the trigger particle 
efficiency cancel out in the normalization. 

Tracks that are spatially near each other can be com- 
bined into a single reconstructed track due to merged 
space points of ionization in the STAR TPC. This track 
merging results in a pair inefficiency at A77 ~ and small, 
but finite A</> whose value depends on the magnetic field 
polarity, charge combination and the pt's of the trigger 
and associated particles [Hj]. The track merging effect 
is most significant in central collisions where the TPC 
hit occupancy is high. The track merging effect in our 
centrality range of 20-60% is negligible. 



A. Event-Plane Reconstruction 

We use the second Fourier harmonic in azimuthal an- 
gle to determine the event-plane angle iJjep [3~lj . which 
is not identical to the real reaction-plane angle due to 
finite resolution. The particles used to determine the 
event plane are below px — 2 GeV/c. To avoid self- 
correlations, particles from the px bin that is used in 
the correlation analysis are excluded from event-plane 
reconstruction. For example, for the associated parti- 
cle pr bin of 1.0 < < 1.5 GeV/c, the particles 
used to calculate the event plane are from px ranges of 
0.15 <p T < 1.0 GeV/c plus 1.5 < p T < 2.0 GeV/c. We 
use the pr-weight method [3l| which gives better event- 
plane resolution due to the stronger anisotropy at larger 
Pt- The slight non- uniform efficiency and acceptance in 
azimuthal angle were corrected as mentioned previously 
in the event-plane reconstruction. Figure [T] shows exam- 
ples of the constructed event plane azimuthal angle distri- 
butions. As seen from the figure, the constructed event 
plane ipEP distribution is approximately uniform. We 
weight the events by the inverse of the event-plane angle 
distributions in Fig. [T] in our correlation analysis. How- 
ever, we found negligible difference in our results with 
and without this event-plane weighting. 
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FIG. 1: (Color online) Constructed event plane azimuthal an- 
gle (V'ep) distributions by the modified reaction- plane (mrp) 
method (points) and the traditional reaction-plane method 
(histogram). The particles used for constructing the event 
plane shown in this figure are from 0.15 < pr < 1 GeV/c or 
1.5 < pr < 2 GeV/c, to be used for correlation analysis for 
the associated particle p£' bin of 1 < < 1.5 GeV/c. 



Non-flow correlations, such as di-jets, can influence the 
determination of the event plane. To reduce this effect, 
we exclude particles within pseudorapidity difference of 
1*7 — f?icad| < 0.5 from the highest pr particle (leading 
particle 77i ca( j) in the event from event-plane reconstruc- 
tion. This method is called the modified reaction-plane 
(mrp) method [32j |. Figure Q] shows in the data points 
a distribution of the constructed t/'ep from this method. 
However, we also use the difference in azimuthal corre- 
lations between analyzed at small and large \Arj\ to ex- 
tract the near-side jet-like component. The MRP method, 
which excludes particles within |A?j| < 0.5 of the lead- 
ing particle in the event, would have different systematic 
biases on the A</> correlations at small and large \Ar]\. 
Thus, for the jet-like component extraction we use the 
traditional reaction-plane method, where no specific cut 
relative to the leading particle is applied in event-plane 
reconstruction. Figure [1] also shows the ^ep distribu- 
tion constructed by the traditional method. We have 
checked the correlation between the event plane angles 
constructed from the traditional method and the MRP 
method, and found they are well correlated. 

We divide our data into six equal-size slices of trig- 
ger particle azimuthal angle relative to the event plane, 
4> s , and analyze azimuthal correlations separately in each 
slice. Figure [2] shows a schematic view, with the slices 
numerically labeled 1 to 6 corresponding to tft s = \<p t — 
ip EP \ = 0°-15°, 15°-30°, 30°-45°, 45°-60°, 60°-75°, and 
75°-90°. We form azimuthal correlations with trigger 
particles in each slice separately. Figure [3] shows, as ex- 
amples, the raw azimuthal correlations in 20-60% Au+Au 
collisions for six slices in 4> s for trigger and associated par- 
ticle pt ranges of 3 < pip < 4 GeV/c and 1 < < 
2 GeV/c (upper panel), and 4 < pip < 6 GeV/c and 
2 < pip^ < 4 GeV/c (lower panel), respectively. All raw 



correlation functions are presented in Figs. [221 123 [Ml 
and [25] in Appendix as a function of trigger pip , asso- 
ciated pi^ , and 4> s . 

Out-of-plane 

A 




Out-of-plane 

FIG. 2: Sketch of six slices in trigger particle azimuthal angle 
relative to the event plane, cp a = \<pt — tpEp\- 



B. Anisotropic Flow Background 

The correlation structure sits atop a large background. 
The background has a flow modulation induced by the 
anisotropies of the trigger particle and the background 
particles with respect to the reaction plane 33]. The 
flow correlated background is given by [34| 



dN 
~dA4> 



= B 



1 + 2v. 



(<*)„,(*,#) 



cos(2A0) + 2vf'vf H) cos(4A0) 

(1) 

where B is the background normalization. In Eq. (JlJ 
and are the associated particle's second and fourth 



harmonics, and arid v^ , " A> are the average harmon- 

ic R) 

ics of the trigger particles, v\' = (cos [2 ((fit — ipBP)]) 



(R) 



and t)^'"' = (cos [4 (4> t — iPbp)])^ R \ respectively. The 
superscript l (R)' indicates that the averages are taken 
within the <f> t region of a slice of width 2c at <p s : 4> s — c < 
\4>t — V'epI < 4>s + c (where c = 7r/24 in our analysis). 
Note that we have used <j) s here and in Eq. §5§ to be 
the center value of a \<f>t — t/'epI bin, while elsewhere we 
simply use <j) s — \4>t ~ V'epI to indicate a narrow bin in 
\4>t — "0EP I ■ For the n th harmonic anisotropy we have [34| , 
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FIG. 3: (Color online) Raw dihadron A<ji> correlations with trigger particles in six slices of azimuthal angle relative to the event 
plane, <f) s = \cj>t — i^ep\- The data are from minimum-bias 20-60% Au+Au collisions. The trigger and associated particle pr 
ranges are 3 < p^) < 4 GeV/c and 1 < p^ ] < 2 GeV/c (upper panel), and 4 < p%? < 6 GeV/c and 2 < pft < 4 GeV/c (lower 
panel), respectively. Note the lower panels correspond to the kinematic range used in Ref. |{|. Both the trigger and associated 
particles are restricted to within \r)\ < 1. Statistical errors are smaller than the symbol size. The curves are flow modulated 
background calculated with four particle i>2{4} and away-side two-particle V2{2,AS} (dashed curves), and the average V2 from 
the two methods (solid curve). 



Here T k is a short-hand notation for 



. sin(fcc) . 
cos(fc0 s ) — £ — -(cos(fcAVO) , 



kc 



(3) 



(cos(kAtp)) = (cosfc( , 0Ep — ipup)) is the event-plane res- 
olution with respect to the k th harmonic, and S is Kro- 
necker's delta. Since the correlation signal we are study- 
ing is on the order of a few percent of the background, 
we need to keep the flow correction in Eq. @ up to the 

0.1%. Keeping terms for n£ up to and 



order V2V4 

for V4 up to i>2, we have 

T 2 + (1 + T^v 



and 



(t,R) 



,(t,R) 



(*) 
2 



(T 2 



1 



2T 4 v\ 



(t) 



(4) 



Ti + (T 2 



(1 



1 + 2T 2 v! 2 t) + 2T A v, 



,(*) 



Ti + (T 2 + T 6 )v. 



(t) 



1 + 2Tou 



(t) 



(5) 



2^2 



The final flow correction is given by Eqs. ([T]), ((4|), and 

The event-plane resolutions, (cos(fcA'!/>YWfc = 2,4,6), 
are obtained from the sub-event method [3l]. The event 
is randomly divided into two sub-events a and b with 
equal multiplicities. The sub-events, excluding the as- 
sociated particle px region, are analyzed to yield event- 
plane angles which, ideally, should be identical. The dif- 
ference between the obtained event-plane angles, ipa — ^b, 
gives the uncertainty in the event-plane determination of 
the sub-events (3l| 

(cOs(fcAVO)sub-cvcnt = V (cOsfc(V>a ~ Ipb)) ■ (6) 



The event-plane resolution of the full event can be ap- 
proximated by [3lT | 



(cos(fcA?/0) « \/2(cos(A;AV0) 



sub— event 



(7) 



in the limit of small event-plane resolution. The factor 
\[2 comes in because the multiplicities of the sub-events 
are smaller than the full event multiplicity by a factor 
of 2. We use the approximate form of Eq. (0) to assess 
systematic uncertainties in the event-plane resolutions by 
different ways of dividing the event into sub-events (see 

sec. mm). 

The precise form of the event-plane resolution of the 
full event is given by [3l| 



(cos(fcAV')) 
where 



X k(N) = v k V2N 



(Pt) 



(8) 



(9) 



depends on the harmonic anisotropy magnitude v k and 
the number of particles TV used in event-plane reconstruc- 
tion. The pt enters into Eq. (j9]) because we weighted each 
particle by its px in constructing the event plane. In data 
analysis we solve for the sub-event Xk(N/2) by Eq. (JSJ 
and the known event-plane resolution of the sub-events 
from Eq. ^ employing an iterative procedure [3lT ] . From 
Eq. © we obtain the full event Xk(N) = V2~Xk(N/2). 
We then use Eq. ([S]) to determine the event-plane resolu- 
tion of the full event 3lJ. The event-plane resolutions are 
listed in Table Q] The resolutions depend on the pr bin 
because particles in a given px bin (to be used for corre- 
lation analysis) are excluded from the event-plane recon- 
struction to avoid self-correlations as aforementioned. 
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One would naively expect that the event-plane res- 
olution should be different for different trigger particle 
orientations from the event plane because the influence 
of di-jets on the event-plane determination should vary: 
a di-jet aligned with the reaction plane enhances the 
event-plane reconstruction resulting in a better resolu- 
tion, whereas a di-jet perpendicular to the reaction plane 
reduces the accuracy of the constructed event plane re- 
sulting in a poorer resolution. However, this is a post 
effect due to the selection of the trigger particle angle rel- 
ative to the event plane. The resolutions used in Eq. (fTJ), 
on the other hand, are those of all triggered events before 
any selection of the trigger particle orientation is made. 
We have also verified this with Monte Carlo toy model 
simulations. 

Since only triggered events enter into our correlation 
measurements, the event-plane resolutions are measured 
using only these events. The event-plane resolutions from 
inclusive events (minimum-bias events within the given 
centrality bin) are found to be within a couple of percent 
of that from the triggered events (see systematic uncer- 
tainty discussion in Sec. IIII Bj) . 

We analyzed the elliptic flow in each of the px bins used 
in our correlation analysis. The obtained elliptic flow 
parameters are tabulated in Table U together with their 
systematic uncertainties. The analysis of the elliptic flow 
and the assessment of its systematic uncertainty are both 
described in Sec. IIII Al We used these v 2 parameters for 
background subtraction. The calculated magnitudes of 
the elliptic flow modulation, 2v < £'' 't^' j are listed in Ta- 
ble fll] together with their systematic uncertainties. The 
calculated background curves are superimposed in Fig. EI 

As mentioned previously, our trigger particle px ranges 
are 3 < < 4 GeV/c and 4 < pf < 6 GeV/c. In 
elementary p+p and ci+Au collisions, the particles in 
these pt ranges originate mainly from hard-scatterings 
and jets. In relativistic heavy ion colllisions, however, 
a large baryon to meson ratio was observed in the pt 
region around 3 GeV/c [H, [36j|. The reason for the 
large ratio and the sources of those high px particles are 
still under debate. The coalescence and recombination 
models [37i - [39| can elegantly explain the large baryon to 
meson ratio from a thermal bath of constituent quarks. 
On the other hand, the jet-like correlations at small an- 
gles relative to trigger particles of p}p > 3 GeV/c, with 
the long range ridge correlation removed, are measured 
to be invariant from p+p, d+Au, peripheral to central 
Au+Au collisions [4l|, and independent of the reaction 
plane direction in Au+Au collisions as will be shown in 
this work. These experimental evidences strongly sug- 
gest that those pj) > 3 GeV/c particles are mostly of jet 
origin in Au+Au collisions just as in p+p and d+Au col- 
lisions. It is possible that recombination may still be at 
work in our trigger particle px ranges, in such that the 
parton(s) prior to recombination have already imprint 
angular correlations related to hard-scatterings (ioj . 

Different sources, such as the recombination (37l - [39| 



and jet fragmentation discussed above, will likely give dif- 
ferent anisotropics to those high px particles. However, 
the anisotropy of the trigger particles to be used in back- 
ground subtraction in Eqs. ([!]), ©, Q, and §5§ should 
still be the experimentally measured net anisotropy [42j , 
as we have done in this work, irrespective of the different 
origins. 

In Eq. (fTJ) we have neglected the odd harmonic terms, 
such as 2v[ a ^v^' R ^ cos(A^) and 2v^v^' R ^ cos(3A0). 
Due to symmetry at mid-rapidity, the averages of the odd 
harmonic coefficients v\, Vq, and etc. vanish. However, 
their fluctuations would yield non-vanishing averages of 
the products of v^v^ and Ug vi , thereby contributing 
to the background in the dihadron correlations. If one as- 
sumes that the amplitude of the v\ (directed flow) fluctu- 
ations is of the same order of magnitude as the maximum 
v i in our pseudorapidity range (which was measured to 
be small [ijiyjthen the v\ fluctuation contribution can 
be neglected [44], |4f| . In the present work we neglect any 
direct flow fluctuation effect in our background subtrac- 
tion. Note that the possible effect of statistical global 
momentum conservation can generate a negative dipole 
which has the same shape as the v\ fluctuation effect. 
However, the statistical momentum conservation effect 
is not from v\ fluctuations, but part of the correlation 
signal, as same as momentum conservation by any other 
mechanisms, such as dijet production. 

Recently, it was argued that the initial fluctuations 
in the overlap geometry (spatial distribution of partic- 
ipating nucleons) may give rise to v% (triangular flow) 
fluctuations [46|, [47], . It was found from the Monte 
Carlo Glauber model that the triangularity due to geom- 
etry fluctuations can be comparable to the magnitude of 
the eccentricity which is connected to the elliptic flow. 
It is thus possible that large triangular flow fluctuations 
can arise which would give triangular peaks in the flow 
background [H, [47|, [49j] . This appears to be the case in 
the AMPT (A Multi-Phase Transport) model and the 
UrQMD (Ultrarelativistic Quantum Molecular Dynam- 
ics) model studied in Ref. [47[ and [HI, respectively. 

Since the orientation of the triangular overlap shape 
due to fluctuations is random relative to the event-plane 
direction, determined by the elliptic anisotropy, the ef- 
fect of any triangular flow is independent of the reac- 
tion plane. In other words, the triangular flow back- 
ground would be proportional to 2v { 3 a) v^ R) cos(3A0) = 

2«3 t>3 cos(3A0) independent of <p s . We may estimate 
the effect of triangular flow fluctuations in our correla- 
tion measurements. The AMPT and UrQMD models in- 
dicate that in the 20-60% centrality range the triangular 
flow fluctuation effect is about 10% of the elliptic flow for 
our trigger and associated px bins, vf/uf ~ 0.1 [47H49| |. 
Experimental data on inclusive two-particle correlations 
at px > 2 GeV/c indicate a ratio of the harmonic coeffi- 
cients also of the magnitude v\jv\ « 0.1 within 20-60% 
centrality 50]. This suggests that the measured third 
harmonic term in the inclusive two-particle correlations 
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TABLE I: Elliptic flow and event-plane resolutions as a function of pr in 20-60% minimum-bias Au+Au collisions. The 
resolutions depend on the pr bin because particles in a given pr bin are excluded from the event-plane reconstruction to avoid 
self-correlations. The errors on V2 are systematic uncertainties. Systematic uncertainties on the resolutions are negligible. 



p T (GeV/c) 


V2 


(cos(2AV>)) 


(cos(4A</>)} 


(cob(6A^)> 


0.15 - 0.5 


0.033 ± 0.001 


0.673 


0.324 


0.127 


0.5 - 1 


0.079 ± 0.003 


0.596 


0.247 


0.082 


1 - 1.5 


0.124 ± 0.004 


0.637 


0.286 


0.104 


1.5 - 2 


0.155 ± 0.006 


0.676 


0.328 


0.129 


2 - 3 


0.176 ± 0.010 


0.704 


0.360 


0.150 


3 - 4 


0.184 ± 0.022 








4 - 6 


0.157 ± 0.024 









(")„,(*,«) 



as a function of p)j?' (in rows) and <j> s = 



TABLE II: The elliptic flow modulation in the correlation background, 2v. 

\4>t — i/'epI (in columns) in minimum- bias 20-60% Au+Au collisions. Both trigger particle pr ranges of 3 < pP < 4 GeV/c 
and 4 < p^' < 6 GeV/c are listed. Quoted errors are systematic uncertainties. Note the significantly smaller systematic 
uncertainties out-of-plane than in-plane. 



p^ (GeV/c) 



O-vr/12 



tt/12 - 71-/6 



tt/6 - tt/4 



7T/4- tt/3 



tt/3 - 5tt/12 



5tt/12 - tt/2 



0.15 - 0.5 
0.5 - 1 
1 - 1.5 
1.5 - 2 
2 - 3 

0.15 - 0.5 
0.5 - 1 
1 - 1.5 
1.5 - 2 
2 - 3 



0.0480 ± 0.0019 
0.1045 ± 0.0056 
0.1726 ±0.0073 
0.2240 ± 0.0108 
0.2605 ± 0.0182 

0.0473 ± 0.0020 
0.1024 ± 0.0057 
0.1697 ±0.0075 
0.2208 ± 0.0109 
0.2572 ± 0.0183 



0.0381 ± 0.0019 
0.0837 ± 0.0053 
0.1376 ± 0.0071 
0.1778 ±0.0102 
0.2060 ± 0.0163 

0.0371 ± 0.0019 
0.0812 ± 0.0054 
0.1338 ± 0.0073 
0.1733 ±0.0103 
0.2012 ±0.0163 



0.0199 ± 
0.0458 ± 
0.0734 ± 
0.0926 ± 
0.1055 ± 

0.0186 ± 
0.0425 ± 
0.0685 ± 
0.0868 ± 
0.0992 ± 



3 < pf 
0.0016 
0.0045 
0.0062 
0.0083 
0.0118 

4 < pP 
0.0017 
0.0046 
0.0064 
0.0084 
0.0118 



< 4 GeV/c 
-0.0030 ± 0.0010 
-0.0016 ± 0.0028 
-0.0071 ± 0.0041 
-0.0143 ± 0.0046 
-0.0208 ± 0.0045 

< 6 GeV/c 
-0.0043 ± 0.0010 
-0.0050 ± 0.0029 
-0.0121 ± 0.0043 
-0.0203 ± 0.0047 
-0.0273 ± 0.0044 



-0.0243 ± 0.0002 
-0.0457 ± 0.0008 
-0.0821 ±0.0013 
-0.1141 ±0.0003 
-0.1388 ± 0.0043 

-0.0254 ± 0.0002 
-0.0487 ± 0.0009 
-0.0864 ± 0.0015 
-0.1190 ±0.0002 
-0.1439 ± 0.0043 



-0.0372 ± 0.0004 
-0.0723 ± 0.0007 
-0.1273 ±0.0008 
-0.1743 ±0.0032 
-0.2099 ± 0.0094 

-0.0380 ± 0.0004 
-0.0749 ± 0.0006 
-0.1309 ±0.0007 
-0.1782 ±0.0031 
-0.2138 ±0.0094 



at low px may be dominated by triangular flow fluctu- 
ations, just as the second harmonic term dominated by 
elliptic flow. As we will show in Sec. IIVCI the effect of 
the triangular flow of this magnitude is relatively small 
in our dihadron correlation measurements with high px 
trigger particles. Because of the lack of systematic mea- 
surements of 7J3 fluctuations as a function of centrality 
and px , in the present work we neglect the possible con- 
tributions from 7J3 fluctuations in our dihadron correla- 
tion measurements relative to high px trigger particles. 

Glauber model calculations also show that the quad- 
rangularity, pentagonality, and hexagonality due to ge- 
ometry fluctuations equal to the triangularity, all large 
and comparable to the eccentricity. However, it was sug- 
gested that those higher order eccentricities were inef- 
ficient to generate sizeable high-order harmonic flow in 
final state momentum space [48j . Experimental data also 
indicate that the magnitudes and fluctuations of 7J4 and 
vq are small relative to the magnitude of vi [o"lT | . Further- 
more, there is no evidence of a large w|, , or vi contribu- 
tion in two-particle correlation measurements |10l , [l3l[T3 |. 



Thus, it is safe to neglect the higher order anisotropic 
fluctuation terms in the flow background of Eq. ([1} . 



C. zyam Background Normalization 

The flow correlated backgrounds given by Eq. ([1]) 
are shown in Fig. [3] as solid curves. The back- 
ground curves have been normalized assuming that the 
background-subtracted signal has Zero Yield At Mini- 
mum (zyam) [l(| [52|. To obtain the ZYAM normalization 
factor, we fold the raw correlation function to within the 
range of < Acf> < ir because of the symmetry of the 
correlation function. We take the ratio of the folded raw 
correlation to the background curve of Eq. fl} , where B 
is set to unity before taking the ratio. We obtain a con- 
tinuous range of the size of 7r/6 where the average ratio 
is the smallest. This smallest average ratio is the nor- 
malization factor B to be used in the flow background of 
Eq. |T]) , which is then subtracted from the raw correlation 
function to obtain the final correlation signal. 
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The background levels can be different for the different 
4> 8 slices because of the net effect of the variations in 
jet-quenching with <j) s and the centrality cuts in total 
charged particle multiplicity in the TPC within < 0.5. 
Thus, in our correlation analysis, the background level B 
is treated independently in individual <p s slices. 

Table IIIII lists the obtained background level B as a 
function of 4> s and p^ in 20-60% Au+Au collisions. Re- 
sults from both trigger particle pr ranges of 3 < pj) < 
4 GeV/c and 4 < p$ < 6 GeV/c are listed. The back- 
ground levels listed are not only for the correlation func- 
tions with the | A?/ 1 < 2 region within our acceptance, but 
also for those in the large Arj region of |A?j| > 0.7. The 
latter is used for the ridge studies (see Sec. IIVC[) . One 
notices that the background level for the lower trigger 
particle Pj) range is slightly larger. This is due to the 
fact that relatively more events contain multiple jets with 
the lower trigger particle pj) and those events are used 
multiple times in our di-hadron correlation analysis [l3j . 

It is worthwhile to emphasize here that our quantita- 
tive results depend on the assumption of the ZYAM back- 
ground normalization. However, as will be discussed in 
Sec. IIIIE1 our qualitative conclusions are independent of 
the zyam normalization. 



III. SYSTEMATIC UNCERTAINTIES 

Background subtraction is the major source of sys- 
tematic uncertainty in our results. The background, 
as given by Eq. ([T]), has three important ingredients: 
the anisotropic flow measurements v 2 and Vi, the event- 
plane resolutions, and the background magnitude B. We 
discuss these systematic uncertainties in Sections IIII A[ 
IIIID1 respectively. They have effects on dihadron cor- 
relation functions presented in Sec. IIV Al and away-side 
correlation widths and magnitudes presented in IIVBI 

We also report results on near-side jet-like and ridge 
correlations in Sec. IIV CI Uncertainties in v 2 and the 
ZYAM background normalization contribute to the un- 
certainties in the ridge correlation results. They do not 
affect the jet-like correlation results in which they can- 
cel. Additional systematic uncertainties arise from the 
assumption of a uniform ridge in Arj, which affects both 
the ridge and jet-like results. These additional system- 
atic uncertainties are discussed in Sec. IIII Fl 

A. Systematic Uncertainty due to Anisotropic Flow 

The anisotropic flow (mainly elliptic flow) background 
which is to be subtracted from the dihadron correlation 
is the anisotro py caused by particle correlations to the 
reaction plane |53l [54[ . There are several measurements 
of elliptic flow; many of them are affected to various de- 
grees by non-flow contributions that are caused by par- 



ticle correlations unrelated to the reaction plane, such 
as resonance decays and jet-correlations. One technique, 
called the event-plane method, is to construct the event 
plane from all charged particles except those of interest 
and then calculate v 2 {ep}= (cos 2(0 — V'ep))/ (cos 2 Aip) 
for the particles of interest, where (cos 2 Aip) is the event- 
plane resolution [3lj . This method is affected by non-flow 
contributions in both sets of particles, those of interest 
and those used to construct the event plane. The v 2 {ep} 
already contains flow fluctuation effects which should be 
included in the jet-correlation background. 

Another method, called the two-particle method, is to 
calculate v 2 {2} = \J (cos2A0) using all particle pairs of 
interest [3l[ . This method is affected by non-flow only in 
the interested particles used for correlation studies. This 
flow parameter also contains flow fluctuation effects. 

The third method, called the four-particle method, is 
to obtain ^{4} from the four-particle cumulant [55j. 
This method is less affected by non-flow from particle 
clustering because the non-flow arising from two particle 
correlations is eliminated, and the non-flow from three 
particle correlations docs not contribute. This method is 
subject to non-flow from higher orders (four-particle cor- 
relation and above) but those contributions are a second 
order effect if compared to non-flow from two particle cor- 
relations [55j. The flow fluctuation will give a negative 
contribution to «2{4}. [11] 

The fourth method is to decompose the low px two- 
particle correlation (the so-called untriggered correlation, 
without the requirement of a trigger particle) into a near- 
angle Gaussian, a dipole, and a quadrupole, and infer 
v 2 {2b} from the fitted quadrupole (56|. The method at- 
tempts to geometrically separate the reaction-plane cor- 
related v 2 from other (i.e. non-flow) correlations (small- 
angle correlations and large-angle dipole). However, the 
method assumes a functional form for those non-flow cor- 
relations, and the decomposition at high pt is difficult be- 
cause of limited statistics and significant jet-correlation 
effects. 

The measured v 2 {2} and ^{mrp} are similar and they 
both significantly overestimate elliptic flow due to large 
contributions from non-flow and fluctuations. While 
flow fluctuation effect should be included in our back- 
ground subtraction, non-flow should be excluded. The 
major component of non-flow is the measured small- angle 
two-particle correlation 57]. Non-flow effect at large 
opening angle (i.e. away-side) is significantly smaller. 
The average (cos2A0) computed from the away side 
of untriggered two-particle correlations, therefore, con- 
tains less non-flow than that computed from all particle 
pairs. Note that, in the absence of non-flow, the average 
(cos 2 Acfi) 's computed from the near-side only, the away- 
side only, or the whole A(j> range are identical. We denote 
the away-side ^/(cos 2A0) |A0|>W2 as v 2{2,as}. If the 
non-flow effect in the away-side untriggered two-particle 
correlation is described by a negative dipole, —Acos Acf>, 
then «2{2,AS} is always larger than the true v 2 {2}: 
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TABLE III: Background level B as a function of p^ (in rows) and (f> s = \<f> t — ipEp\ (in columns) in minimum-bias 20-60% 
Au+Au collisions. Both trigger particle pr ranges of 3 < pp 1 < 4 GeV/c and 4 < p^' < 6 GeV/c are listed. The trigger and 
associated particles are within \r)\ < 1. Backgrounds are tabulated for the entire \Ar/\ < 2 range of our acceptance as well as 
for large Arj cut of \Arj\ > 0.7. The first error is statistical. The second error is the quadratic sum of the zyam systematic 
uncertainty and the one-sided systematic uncertainty due to background deviation from zyam. The former is assessed by varying 
the A<f> normalization range. The latter is assessed by comparing our zyam background to those obtained from asymmetric 
correlations of the separate positive and negative (f>t — fp~EP regions. 



P T 



£° (GeV/c) 



0- tt/12 



tt/12 - tt/6 



tt/6 - tt/4 



tt/4 - tt/3 



tt/3 - 5tt/12 



5tt/12 - tt/2 



0.15 - 0.5 
0.5 - 1 
1 - 1.5 
1.5 - 2 
2 - 3 

0.15 - 0.5 
0.5 - 1 
1 - 1.5 
1.5 - 2 

2 - 3 

3 - 4 

0.15 - 0.5 
0.5 - 1 
1 - 1.5 
1.5 - 2 
2 - 3 

0.15 - 0.5 
0.5 - 1 
1 - 1.5 
1.5 - 2 

2 - 3 

3 - 4 



3 < < 4 GeV/c 

47.23 ± 0.02 
22.92 ± O.Oll ,;™ 23.31 ± 0.011°: 



47.26 ± 0.01 
22.37 ±0.01 



+0.01 
-0.06 
+0.03 
-0.09 

5.986 ± 0.0051°, ° ll 
1.662 ± 0.0021°;^ 
0.642 ± 0.002l l ;° )3 1 



46.47 ± 0.04 
22.07 ±0.02 
5.909 ±0.012 
1.644 ± 0.006 



+ 0.05 

0.11 
+0.03 
0.06 
+0.010 
0.048 
+0.006 
0.041 

0.604 ± 0.0041°;°°;? 

0.058 ± o.ooil^^a 



19.31 ± O.OII0.05 
9.14 ± 0.0ll°;°i 
2.437 ± 0.003lo.°°g 
0.674 ± 0.002l°;°°4 
0.258 ± 0.0011o.to4 



18.92 ± 0.02l°;°J 
8.98 ± 0.0ll°,;°j 4 
2.406 ± 0.008l^?o 
0.666 ± 0.0041°;°^ 
0.239 ± 0.0021° °°^ 
0.023 ± 0.00ll°'o°i 



47.26 ± 0.021°;°^ 
22.58 ±0.0ll°,;?5 
6.039 ± 0.0051^-^8 
1.672 ± 0.0021° °°^ 
0.651 ± 0.002lg;°^ 



46.56 ± 0.041° !° 
22.23 ± 0.021°,;?,! 
5.958 ± 0.012l° °°* 
1.651 ± 0.006lo °°5 
0.612 ± 0.0041°;°™ 
0.059 ± O.OOll ,; , ^ 



19.32 ± O.Oll^EJg 

9.24±O.OllHo 
2.461 ± 0.003lg;^g 
0.681 ± 0.002l°;°°J 
0.264 ± 0.0011° °^ 



18.94 ± 0.02l° °J 

9.06 ±0.02 
2.422 ± 0.0081°,;°,°,!? 
0.664 ± 0.0041°,;°$ 
0.247 ± 0.003l^;^9 
0.024 ± 0.0011°'°°° 



,+0.05 
-0.12 



47.29 ± 0.021°, % 



+0.06 
0.08 
0.03 
18 

6.130 ± 0.0051^08 6 - 181 ± O-OOStaoS 
1.696 ± 0.003lg;° 4 ° 1.703 ± 0.0031°,;°,°;; 



4 < p^> < 6 GeV/c 

46.77 ± 0.04 
22.41 ± 0.02l^S 23.12 ± 0.03 



46.72 ± 0.041°,;°!? 



5.986 ± 0.0131!!'??!! 6.117 ± 0.013 



0.610 ± 0.004lg;g°2 0.623 ± 0.004l°;° 2 ° 
0.061 ±0.001 



3 < p ( * ] < 4 GeV/c, |Ar?| > 0.7 



19.32 ± 0.0ll°,;°j 
9.36 ± O.Oll ,;? , 
2.497 ± 0.0031°,°,°,? 



19.01 ±0.03 



+0.02 



18.97 ±0.03 



+0.03 



9.11 ±0.02l°°° 
2.436 ± 0.0081°, °£l 2.458 ± 0.008lg;^5 



9.37 ± 0.021°°? 



-0.078 

0.678 ± 0.004l°;°J° 

0.001 
0.014 



0.247 ± 0.003Z 



0.681 ± 0.004 
0.253 ± 0.003 



0.008 
-0.039 
0.001 
0.014 



0.025 ± 0.0011°;°°° 0.025 ± 0.00ll°'°°i 



47.16 ±0.021°,;?* 
23.40 ± O.Oll° °o 
6.130 ±0.005+S-^o 



47.09 ± 0.02 
23.59 ± 0.01 
6.198 ±0.006 



0.08 
0.05 
+0.06 
-0.09 
+0.024 



0.662 ± 0.0021° °°^ 0.663 ± 0.0021° °°^ 

,(*) 



1.695 ± 0.003lo;°°g 1.693 ± 0.0031°,;°^ 
0.655 ± 0.002l° °?l 0.657 ± 0.002l°°°i 



0.06 
0.27 
+0.05 
0.19 
0.037 
0.101 
0.023 
0.051 



+0.07 
0.08 
0.13 
0.08 
0.032 
0.059 
0.018 
0.022 



+0.002 A Ofi o i Q nm +0.000 
0.014 <J.UOO =C U.UU1_ 0>016 



19.25 ±0.0ll^g 
9.51 ± 0.011°;°° 
2.514 ± 0.0031°,;°,^ 

r0.004 

-0.017 

0.269 ± 0.0011°,;°,?°, 0.268 ± O.OOll ,; ,?,! 
4 < pP < 6 GeV/c, |At?| > 0.7 



46.67 ±0.05 

23.07 ±0.03 
6.077 ±0.014 
1.674 ±0.007 
0.616 ± 0.0041°;°^ 
0.061 ± 0.00llg;^J 



19.27 ±0.011° °^ 
9.55 ± 0.011°;™ 

-0.009 
-0.026 
-0.008 
-0.015 



46.75 ± 0.05l°;°* 
23.41 ±0.03l°;ig 
6.171 ±0.014l°°5* 
1.711 ±0.007l°;°°g 
0.614 ± 0.0051°;°?* 
0.061 ± 0.00ll°,;°,°, 4 



2.494 ± 0.003 
0.686 ± 0.002 



19.18 ±0.01 
9.60 ±0.01 
2.517 ±0.004 
0.684 ± 0.002 



+0.01 
-0.02 
+0.03 
-0.03 
+0.012 
0.016 
+0.006 
0.004 



0.265 ± O.OOll ,; ,?^ 0.264 ± COOll ,' ^ 



O.OO.'i 



o.oos 



18.95 ± 0.03 

9.39 ± 0.02 
2.458 ± 0.009 
0.673 ± 0.005 
0.245 ± 0.003 
0.023 ± 0.00ll°;°° 7 



+0.01 
0.03 
+0.01 
0.08 
h0.005 
-0.034 
0.004 
-0.008 
0.003 
0.007 



18.96 ±0.03] 



+0.00 



9.50 ± 0.02 
2.497 ±0.010^ 
0.686 ± 0.005^ 
0.247 ±0.003 



0.024 ± 0.00ll°'°° 2 



+0.003 



^{2,AS} = 



lTi2 (1 - -4 cos + 2v\ . truc {2} cos 2 A<j>) cos 2A(f>dA<p 

ill* ^ ~ A COS A< ^ + 2 «2{2}truc COS 2A0) dA<j) 



J 2,truc 



{2}± 



2A 



1 



7T ± 2 A V 3 



4truc{2} ) > 4true{2}- 



If the away-side non-flow effect is described by a Gaussian, Aexp 
Gaussian is relatively narrow: 



(A0-x) 2 
2cr 2 



, then the same holds as long as the 



«2{2,AS } 



r 37r/2 /, . 


(A0-7T) 2 " 

2<T 2 


± 2tj2 truc {2} cos 2A<p) cos 2A<pdA(f> 


$1* (l + ^^xp 


(A4>- 
2a 


-it) 2 ' 
5 


+ 27j2{2} trU oCos2A(/)^ dA4> 



«2,true{2} ± 



2irAa 



it ± \J2jtAo 



(f^ - vl truc {2}) > v\ truc {2}. 
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In principle non-flow contribution can be negative; 
however, untriggered two-particle correlations show no 
evidence of a negative non-flow [57| . Therefore we take 
«2{2,AS} as the upper limit of the systematic uncertainty 
of the elliptic flow parameter in our background sub- 
traction. Note that «2{2,AS} automatically includes the 
positive flow fluctuation effect. We compute «2{2,AS}= 
J (cos2A0)|A0|>7r/2 using all pairs of particles in each 
given p^ bin with opening angle |A</>| > ir/2. We use 
minimum-bias events within our 20-60% centrality bin 
(i.e. inclusive events) to compute v 2 {2,AS}. Thev2{2,AS} 
computed from triggered events only (with at least one 
trigger particle) is nearly identical to that from inclusive 
events. We also compute «2{2,AS} excluding particles 
within | Aij\ < 0.5 of the highest pr particle in the event 
to reduce the non-flow effect. We find the «2{2,AS} ob- 
tained in this manner is not significantly different from 
that calculated without excluding those particles, as com- 
pared to the uncertainty on v 2 (see below). 

We have also calculated «2{2,as} separately for the 
four 10%-size centrality bins. We took the average 
w|{2,AS} weighted by the number of trigger particles in 
addition to the number of associated particle pairs in 
each centrality bin. The obtained average u|{2,AS} is 
generally smaller than but close to that obtained from a 
single 20-60% centrality bin. We used the latter to be on 
the conservative side. 

For high pt trigger particles, the statistical error on 
«2{2,AS} becomes large. We therefore use V2{ep}, in- 
stead of W2{2,AS}, as our upper systematic bound for 
trigger particle v%. The ^{ep} is measured in 10%-size 
centrality bins. We used the average ^{ep} weighted by 
the number of trigger particles. 

The measured «2{4} likely underestimates elliptic flow 
because the flow fluctuation effect reduces i>2{4} [55j. We 
note that «2{4} may still contain some non-flow effects. 
However, the agreement between v 2 {4} and the elliptic 
flow measurement using the Lee- Yang-Zero method sug- 
gests that such non-flow effects are small ;58] . We there- 
fore use i>2{4} as our lower bound of v 2 systematic uncer- 
tainty, as same as in Refs. We use a parameteri- 
zation of the measured «2{4} as a function of pt in each 
10%-size centrality bin. We first calculate the average 
v 2 {4} weighted by the associated particle multiplicity at 
each px?' within a given p^ bin for each centrality bin. 

We then take the average i>2{4} for each pf' bin from the 
four 10%-size centrality bins weighted by the number of 
trigger particles within 3 < p^) < 4 GeV/c. We used 

the same t>2{4} for 4 < p^p < 6 GeV/c because the aver- 
age i>2{4} obtained using the number of trigger particles 

within 4 < pj? < 6 GeV/c as the weight is practically 
the same. 

As the default u 2 , we use the average between w 2 {2,as} 
and t>2{4} for associated particles. We use the range 
bracketed by V2{2,AS} and i>2{4} as our systematic un- 
certainty on v 2 . For trigger particles we use the average 
between d 2 {ep} and v 2 {4} as the default «2 and the range 



bracketed by ^{ep} and v 2 {4} as its systematic uncer- 
tainty. Table U lists the default v 2 values together with 
systematic uncertainties for different pr bins in 20-60% 
Au+Au collisions. 

We parameterized the V4 measurement [32[ as V4 = 
l.l5v 2 and used this parameterization for both trigger 
and associated particles in our flow correction [2l|. The 
uncertainties in v 2 are propagated to W4. Note that the 
V4 fluctuation effects that is related to the second har- 
monic event plane, which should be included in our flow 
background, are already included in the 114 measurement 
which was carried out with respect to the second har- 
monic event plane [H, |5l|. Fluctuations in v 4 related 
to the fourth harmonic event plane could be potentially 
not small [51j and are not included in the available mea- 
surement of V4, however, these fluctuation effects come 
into our two-particle correlation background as v\ (not 
through the cross-term of v 2 V4) and are therefore negli- 
gible. 

The flow backgrounds are shown by the solid curves in 
Fig. 131 The systematic uncertainties due to anisotropic 
flow parameters are shown by the dashed curves. The 
normalization of each background curve is adjusted by 
ZYAM to match the raw correlation function such that 
the background-subtracted correlation is zero at mini- 
mum (see Sec. Ill C|) . As seen from the figures, the dashed 
curves are not symmetric about the solid curve. This is 
mainly due to the zyam normalization as the normaliza- 
tion region is around A<p ss ±1, not at ±7r/2. 

The coefficient v^V^ in Eq. ([1]) determines the size 
of the modulation in the flow background. These coeffi- 
cients are tabulated in Table HT1 For in-plane trigger par- 
ticles, i?2 is positive as given by Eq. @ or (Q} . The cor- 
related elliptic flow uncertainties in and gives 
a large uncertainty in v 2 a ^v 2 ' R \ For out-of-plane trig- 
ger particles, however, v 2 ' R ^ is negative. The correlated 
uncertainty in and v 2 ' R ^ tends to cancel each other, 

resulting in a small uncertainty in t)j tij ■ This is ap- 
parent in the systematic uncertaities listed in Table [TT1 It 
is shown in the systematic uncertainty background curves 
in Fig. [21 where the uncertainty for in-plane correlations 
is large, while for out-of-plane correlations it is small. 

B. Systematic Uncertainty due to Event-Plane 
Resolution 

The event-plane resolutions enter into the flow back- 
ground modulation together with the anisotropic flow pa- 
rameters, via w n (cos(A;A'0)). Terms with k = n are not 
affected by uncertainties in the event-plane resolutions, 
because w n (cos(nA-0)) are the measured anisotropic flow 
parameters. The event-plane resolutions (cos(kAip)) of 
different A;'s are likely correlated, hence the uncertainty 
in v n (cos(kAip)) for k 7^ n due to uncertainties in the res- 
olutions may be greatly reduced. To be conservative, we 
assume the uncertainties in the event-plane resolutions 
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to be uncorrelated in our estimation of their effects on 
our correlation results. 

The systematic uncertainty of the event-plane resolu- 
tion was determined by repeating the sub-event method, 
but splitting the particles by charge instead of by ran- 
dom determination, as done in the default case. They 
are also assessed by comparing the event-plane resolution 
from triggered events only (default) to inclusive events, 
and by applying a weighting of the number of trigger 
particles (default) and not applying this weighting. In 
addition, differences in event-plane resolutions are as- 
sessed with (default) and without event-plane flattening 
by weighting of the inverse of (^-dependent efficiencies. 
The event-plane resolution uncertainties thus estimated 
are typically less than 1% for (cos(2A^>)}, and less than 
2-3% for (cos(4AV0) and (cos(6AV>))- 

The effects of the estimated event-plane resolution un- 
certainties on the final background-subtracted correla- 
tion functions are significantly smaller than those caused 
by the uncertainties on anisotropic flow, and are therefore 
neglected. 

C. Effect of Finite Centrality Bin Width 

For the data reported in this paper, the entire 20-60% 
Au+Au centrality range is treated as a single central- 
ity bin in which the event-plane resolutions and ellip- 
tic flow are obtained and the azimuthal correlation is 
analyzed. Alternatively, the analysis was repeated in 
each of the four 10%-size centrality bins using the cor- 
responding event-plane resolutions and the elliptic flow 
measurements. Those correlation results are added to- 
gether, weighted by the number of trigger particles in 
each centrality bin. The recombined results are consis- 
tent with using a single 20-60% centrality bin, well within 
the systematic uncertainties due to those in flow sub- 
traction and ZYAM normalization. This is because the 
measured elliptic flow V2 is fairly constant over the entire 
20-60% centrality range, so that (v^v^) (v^){v^ a) ). 
The event-plane resolutions vary with centrality mainly 
due to the multiplicity change. However, the event-plane 
resolutions enter into the flow background of Eq. (JT|) lin- 
early, and because the high px trigger particle multiplic- 
ity scales almost linearly with the total multiplicity, the 
effect of the centrality-varying event-plane resolution is 
minimal in the flow correction calculated from the single 
20-60% centrality bin or summed from multiple narrower 
centrality bins. 

D. Systematic Uncertainty due to ZYAM 
Background Normalization 

Naively one would expect the background level B in 
Eq. ([T]) to be the same for all <j> s slices because the un- 
derlying background should not depend on the signal (or 
orientation of the trigger particle). However, there could 



be biases in the event samples with trigger particles at 
different <fi s such that they contain slightly different un- 
derlying background multiplicities due to the possible dif- 
ference in jet-like correlated multiplicities at different <fi s 
and the overall constraints caused by centrality cuts on 
the reference multiplicity. In our analysis we use differ- 
ent B values for different <fi s slices, each independently 
obtained using ZYAM on the correlation function of the 
corresponding slice. 

One source of systematic uncertainty on B is due to the 
limited range in A<fi where the background-subtracted 
correlations appear to have a minimum 'plateau'. This 
part of the systematic uncertainty is assessed by varying 
the size of the normalization range in Acfi between 7r/12 
and 7r/4 (default range is 7r/6), similarly to Ref. [Toj . 

The ZYAM assumption likely gives an upper limit to 
the underlying background level. One could make an 
improved assessment of the background level with more 
stringent requirements, such as using three-particle corre- 
lation ZYAM [2l|. However, the analysis of three-particle 
correlation within a limited <fi s range of the trigger par- 
ticle is difficult. 

In this paper, we assess this part of the systematic un- 
certainty on B by comparing to the ZYAM backgrounds 
obtained separately from correlation functions at posi- 
tive (fit ~ ^ep and negative (fit — V'ep • Those ZYAM back- 
grounds are always lower than our default B from zyam 
of the combined correlation function of positive and neg- 
ative (fit — ipEP- This is because the separately analyzed 
correlation functions are asymmetric about A<fi = and 
Acfi = 7r, and the ZYAM is determined by only one side 
of the correlation function [59l . |60| , whereas in our com- 
bined correlation functions reported here, the two sides of 
the separately analyzed asymmetric correlation functions 
are averaged. We treat the difference between the zyam 
background from this paper and that obtained from the 
asymmetric correlation functions as an additional, one- 
sided systematic uncertainty on B. 

We may also study the background level by fitting the 
ZYAM-background-subtracted correlation functions with 
a combination of Gaussians and a free parameter for an 
offset from zero. Specifically we fit the correlation data 
to three Gaussians (a near-side Gaussian at Acfi = and 
two away-side Gaussians symmetric about Acfi — 7r), and 
four Gaussians (adding a fourth Gaussian at A<fi = ti 
with the same width as the near-side Gaussian). Some 
of the fits yielded unphysical offsets because of the limited 
constraint of the correlation data on the fit model. For 
the other fits, the fitted offsets are comparable to the 
systematic uncertainty obtained from the comparisons 
to the asymmetric correlation functions discussed above. 
The Gaussian fits to the correlation functions without 
the offset will be discussed in Sec. IIV El 

The different sources of systematic uncertainties on B 
are added in quadrature. The total systematic uncer- 
tainty is listed in Table Mil together with the statistical 
uncertainty. We take the quadratic sum of the statistical 
and systematic uncertainties as the total uncertainty for 
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B on our correlation results. 



E. Is the Away-Side Double-Peak an Artifact of 

ZYAM? 

As will be shown in IIV Al the background-subtracted 
correlation functions on the away side are single-peaked 
at Acj) = it for triggered particles in-plane, but double- 
peaked for trigger particles out-of-plane beyond the 
flow systematic uncertainties. Since the subtracted 
background is flow-modulated, the natural question is 
whether the away-side double-peak structure is due to 
an unrealistic systematic uncertainty. To address this 
question, it is worthwhile to note that the flow back- 
ground modulation changes phase when the trigger par- 
ticle moves from in-plane to out-of-plane, as shown in 
Fig. [3l A smaller elliptic flow would make the in-plane 
correlation more peaked at Acj) = and 7r and the out-of- 
plane correlation more dipped at n (hence more double 
peaked on the away side). On the other hand, a larger 
elliptic flow would make the out-of-plane away-side corre- 
lation less double-peaked. One would need a ~15% larger 
v^v^' 1 ^ than in Table [TTl significantly beyond the sys- 
tematic uncertainty from the anisotropy measurements, 
to eliminate the away-side double-peak for the out-of- 
plane 4>s slice. However, this large t4*'^ would result 
in double-peaked away-side correlations for some of the 
other <j) s slices. 

The background magnitude affects the absolute mag- 
nitude of the flow modulation subtracted from the raw 
data in obtaining the correlation signal. Since the back- 
ground normalization is determined by the zyam descrip- 
tion, the question arises whether the away-side double- 
peak for the out-of-plane cp s slices is an artifact of a sig- 
nificantly smaller background level than zyam beyond 
the ZYAM normalization systematic uncertainty. The an- 
swer is negative because the flow background is the low- 
est at Acj> = 7r for out-of-plane trigger particles. Allowing 
a non-zero flow-modulated "pesdestal" into the correla- 
tion signal will exaggerate the double-peak feature, i.e., 
the dip at Acj> = n will be even deeper than the double 
peaks. In other words, if the true background is lower 
than ZYAM, then the away-side correlation functions for 
out-of-plane trigger particles will be more double-peaked. 
Only when the background is larger than ZYAM would the 
dihadron correlation signal become single-peaked; how- 
ever, as a result the signal strength would become nega- 
tive. 

In summary, to eliminate the away-side double-peak, 
one needs either a larger anisotropic flow than measured 
while fixing the background normalization by zyam, or 
a larger background normalization than ZYAM while fix- 
ing the anisotropic flow as measured. To investigate fur- 
ther the interplay between background normalization and 
anisotropic flow and its effect on the dihadron correlation 
signal, we performed a study of free fits to the raw cor- 
relation data, treating the anisotropic flow and the back- 



ground magnitude as free parameters. In order to do so, 
one needs a prescription for the correlation signal func- 
tional form. It has been shown that the sum of a near-side 
Gaussian, a negative dipole, and a quadrupole (reflect- 
ing elliptic flow) can adequately describe the two-particle 
azimuthal correlation at low pt without the requirement 
of a high pt trigger particle [571 ] . Thus, we fit our raw 
correlation data by 



— — = 5(l + 2V r 2 cos2A0 + 2F 4 cos4A0) + 
dAcp 

A ns exp (- ) ~ ^dipoic cos A(f>, (10) 



treating the flow modulations V2 and V4, the near-side 
Gaussian parameters A ns and cr ns , and the negative 
dipole magnitude Adi po ie as free parameters. Figure |31 
(upper panels) shows the fits by Eq. (fTO)) to the raw cor- 
relation functions in six cj) s slices for 3 < Pj) < 4 GeV/c 
and 1 < p^' < 2 GeV/c. The fits are shown in the solid 
curves. The dashed curves show the fitted flow back- 
grounds. The lower panels of Fig. [4] show the correlation 
functions after subtracting the fitted flow backgrounds. 
The fitted near-side Gaussian and the negative dipole are 
depicted individually. 

As seen from the x 2 / NDF written in each upper panel, 
the fits by Eq. (flT))) are generally good. This is also 
true for the other pj) and p^? bins. However, the fit- 
ted flow modulations (written in the lower panels) are 
significantly larger than the measured ones for the out- 
of-plane 4> s slices, much beyond their systematic uncer- 
tainties quoted in Table |IIJ In other words, in order to 
eliminate the away-side double-peak, an anisotropic flow 
that is much larger than experimentally measured is re- 
quired, consistent with our earlier observation. More- 
over, the deviations of the fitted flow modulations from 
the measured ones vary from slice to slice, which should 
not be the case if the measured flow parameters that we 
used were simply in error. Qualitatively the same fea- 
tures are observed for the other pj) and p^ bins. These 
free fit results demonstrate that the fit model of Eq. (JTDJ) 
cannot be the correct functional form to describe the di- 
hadron correlation signal with a high px trigger particle. 

We have also used other single-peaked functional 
forms, e.g., a near-side Gaussian and an away-side Gaus- 
sian, in our fit. Similar conclusions were obtained. The 
away-side double-peak for the out-of-plane trigger parti- 
cles cannot be eliminated without using a flow subtrac- 
tion much larger than experimentally determined, either 
with or without zyam. Thus, we conclude that the away- 
side double-peak structure is not an artifact of the zyam 
flow subtraction procedure used in this analysis. 
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FIG. 4: (Color online) Fit by Eq. (|10p to the raw correlation data in the upper panels of Fig. [3] (Upper panels) The solid 
curves are the fit results and the dashed curves are the fitted flow background. (Lower panels) The correlation functions after 
subtracting the fitted flow background. The fitted same-side Gaussian and negative dipole are depicted individually in the 
dashed and solid curves, respectively. 



F. Systematic Uncertainties on Jet-Like and Ridge 
Correlations 



To obtain the jet-like component we take the differ- 
ence of the correlation functions from \Ar/\ < 0.7 and 
|A?7| > 0.7 (properly weighted by the two-particle Arj 
acceptance). The assumption in this procedure is that 
the ridge is uniform in A77 (after taking into account the 
trivial two-particle A77 acceptance) and is therefore sub- 
tracted away in the difference [ll[ . Measurements at low 
Pt without a trigger particle indicate that the ridge is 
broad but drops with increasing Arj 61|. If this is true 
for trigger particle correlations as studied here, our "jet" 
measurement contains a residual ridge contribution. To 
estimate this effect, we study Arj correlation functions 
for near-side associated particles (|A</>| < 1). An ex- 
ample is shown in Fig. [3] for 3 < pip < 4 GeV/c and 

0.15 < p£° < 3 GeV/c in the 20-60% centrality bin. The 
in-plane direction (0° < <fi < 45°) is used because, as 
will be shown later, the ridge resides mainly in the in- 
plane direction. We compare the ridge contributions to 
the |A?7| < 0.7 region as extrapolated from a constant 
ridge fit and from a linear fit [llj , both done in the large 
Ar] range of \Ar]\ > 0.7. Because of possible edge effects 
in the Arj acceptance, we also limit our fit range within 
0.7 < |A?/| < 1.6. We assign the difference, ±15%, as the 
systematic uncertainty on the jet-like component yield 
due to the assumption of the uniform ridge. 

In this paper, we consider all correlated particles at 
|A?7| > 0.7 to be part of the ridge. The ridge yield we 
report in this paper is defined to be the integral of the 
correlated particle yield over 0.7 < |A?7| < 2.0. Thus, the 
assumption of the ridge shape does not affect the ridge 
yield. 

We have assumed that the jet-like component is con- 
tained within |A?7| < 0.7, and assigned the entire cor- 
related yield in IA77I > 0.7 as ridge. This introduces 



uncertainty in the ridge yield as well as in the jet-like 
yield. Moreover, the fraction of the jet-like component 
that leaks out of the Arj cut is subtracted in obtaining 
the jet-like part, thus the effect of the leakage is doubled 
in the extracted jet-like component. To study this effect, 
we fit the Arj correlation function (such as that shown in 
the upper panel of Fig. [SJ to a Gaussian with centroid 
at A?y = and a constant pedestal (i.e., uniform ridge). 
The Gaussian width is shown in Fig. [5jb) as a function of 

4> s for 1 < p^ < 2 GeV/c and in Fig.[5^c) as a function 

of PjT for integrated (j) s . The Gaussian width does not 

significantly depend on <p s or . We estimate the ef- 
fect of the leakage of the jet-like component to be about 
10% of the jet-like yield, assigned as a single-sided neg- 
ative uncertainty on the ridge yield, and a single-sided 
positive uncertainty, twice as large, on the jet-like yield. 
The physics of the correlation widths will be discussed in 
Sec. lIFEl 

The systematic uncertainty on the jet-like yield due to 
flow uncertainty is small because the large uncertainties 
due to V2 are cancelled, assuming vi is constant over Arj. 
This should be a good assumption because the PHO- 
BOS experiment found that V2 was constant within the 
j] acceptance of the STAR TPC (dropping only towards 
larger [63,l63||. 

Figure [5] illustrates the various systematic uncertain- 
ties on the extracted ridge yield, (i) The systematic un- 
certainties due to flow subtraction are shown by the solid 
curves. The uncertainty is dominant at small <p s ; the 
flow uncertainty at large (j) s is small, (ii) The system- 
atic uncertainty due to background normalization uncer- 
tainty is shown in brackets, as assessed by varying back- 
ground normalization range and by comparing to back- 
ground normalizations of asymmetric correlation func- 
tions at positive and negative <f>t — V'ep separately, (iii) 
There is an additional systematic uncertainty in the ex- 
tracted ridge yield because the jet-like correlation can 
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FIG. 5: (Color online) (a) Raw Ar] correlation of near-side 
associated hadrons (|A</>| < 1) integrated over 0° < (f> a = 

\<jH - itep\ < 45° and 0.15 < p£° < 3 GeV/c, corrected 
by the two-particle Ar] acceptance. The dotted curve is a 
single Gaussian fit and the dot-dashed horizontal line is the 
fit pedestal; the solid lines are linear fits to the regions 0.7 < 
| Arj\ < 1.6 and 0.7 < A77 < 2.0, respectively, and the dashed 
lines are their extrapolations, (b) Gaussian fit a to near-side 
Ar] correlation in 0.15 < p£ ' < 3 GeV/c as a function of 
4> s - (c) Gaussian fit a to near-side Ar] correlation integrated 
over 0° < 4> s < 45° as a function of j4? ' . The data are from 
minimum-bias 20-60% Au+Au collisions. The trigger particle 
Pt range is 3 < p^! < 4 GeV/c. Error bars are statistical. 



be broader than 0.7 in A77 and the jet-like yield beyond 
|Afy| > 0.7 is contained in the extracted ridge yield. This 
part of the systematic uncertainty is shown by the ar- 
rows. The total systematic uncertainties are obtained by 
the quadratic sum of the individual sources and shown 
by the boxes. 
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FIG. 6: (Color online) Illustration of the different systematic 
uncertainties on the ridge yield (defined in Sec. IIV C[) within 
\A<fi\ < 1 and |A?7| > 0.7 as a function of <f) a — \<j> t — "4>ep\- 
The data are from minimum-bias 20-60% Au+Au collisions. 
The trigger and associated particle pr ranges are 3 < pip < 
4 GeV/c and 0.5 < Pj! ' < 1 GeV/c, respectively. Statistical 
errors are smaller than the symbol size. The various system- 
atic uncertainties are due to: (i) flow subtraction (shown by 
the solid curves), (ii) background normalization uncertainty 
(shown in brackets), assessed by varying background normal- 
ization range and by comparing to zyam from asymmetric 
correlations separately for positive and negative (fit — ipEP, and 
(hi) leakage from the jet-like component into the \Ar]\ > 0.7 
region (shown by the arrows). The total systematic uncer- 
tainties are shown by the boxes. 



IV. RESULTS AND DISCUSSIONS 

A. Correlation Functions 

Figure [7] shows the background-subtracted dihadron 
azimuthal correlations in 20-60% Au+Au collisions as 
a function of the trigger particle orientation relative to 
the event plane, <f> s . The thin histograms embracing the 
shaded area indicate the systematic uncertainties due 
to anitropic flow. The horizontal shaded band around 
zero indicates the systematic uncertainties due to ZYAM 
background normalization. The slight modulations of 
the edges of the band are because of the anisotropic 
flow in the combintorial background. For comparison 
the minimum-bias d+Au inclusive dihadron correlation 
(without differentiating with respect to an "event plane" ) 
is superimposed in each panel in Fig. [7] The trigger and 
associated particle pr ranges are 3 < pj) < 4 GeV/c 
and 1 < p^ ' < 2 GeV/c (upper panel), and 4 < p)p < 

6 GeV/c and 2 < p^ < 4 GeV/c (lower panel), respec- 
tively. These kinematic ranges correspond to those for 
the raw correlations shown in Fig. [31 The background- 
subtracted correlations for all trigger and associated par- 
ticle pt ranges are presented in Appendix [5] in Figs. [26l 
and[271 

As seen in Fig. [3 the near-side peaks in Au+Au colli- 
sions are evident for all trigger particle orientations. The 
single-peak shape of the near-side correlation remains rel- 
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FIG. 7: (Color online) Background-subtracted dihadron correlations with trigger particle in six slices of azimuthal angle from 
the event plane, cj>s = \<j>t — ^ep|- The trigger and associated particle pr ranges are 3 < pip < 4 GeV/c and 1 < p^? ' < 2 GeV/c 
(upper panel), and 4 < p!p < 6 GeV/c and 2 < < 4 GeV/c (lower panel), respectively. Note the bottom row corresponds 
to the kinematic range used in Ref. ((J. Both the trigger and associated particles are restricted to be within \r)\ < 1. The data 
points are from minimum-bias 20-60% Au+Au collisions. Systematic uncertainties are shown in the thin histograms embracing 
the shaded area due to flow subtraction and in the horizontal shaded band around zero due to zyam background normalization. 
Statistical errors are smaller than the symbol size. For comparison, the inclusive dihadron correlations from d+Au collisions 
are superimposed as the thick (green) histograms (only statistical errors are depicted). 



atively unchanged from in-plane {<j) s ~ 0) to out-of-plane 
(4> s ~ 7r/2). However, the amplitude of the near-side 
peak decreases with <fi s , becoming similar to that from 
d+Au collisions at large 4> s . Our previous studies have 
shown that the near-side correlation, while not much 
modified at high px, is enhanced in Au+Au collisions 
relative to P+P and d+Au collisions at low to modest 
Pt EH EL EH ■ The present results show that the near- 
side enhancement is mostly present for trigger particles 
oriented in-plane and the modification for trigger parti- 
cles oriented at (f> g ~ ir/2 is minimal for this centrality 
bin. 

Unlike the near side, the away-side correlation struc- 
ture evolves dramatically when trigger particles move 
from in-plane to out-of-plane for the 20-60% centrality 
bin. The away side has a single peak when the trig- 
ger particles are oriented close to the event plane. Only 
when the trigger particle direction is far away from the 
event plane, the double-peak structure emerges on the 
away side. In addition, the away-side modification in- 
creases with increasing associated particle p^ ■ Our pre- 
vious studies showed that the away-side correlation struc- 
ture is significantly modified in central Au+Au collisions, 
and the modification is the largest in the intermediate 
Pt range [h| EH • The present result indicates that the 
away-side modification has a strong dependence on the 
trigger particle direction relative to the event plane. The 
strongest away-side modification is found for trigger par- 
ticles perpendicular to the event plane (see Fig. [7]). The 
results indicate the important role the medium path- 
length traversed by the away-side parton at different 
4> s plays, and should provide useful input to theoretical 
modeling of partonic energy loss in the nuclear medium. 



The lower panel of Fig. [7] shows the high px associated 
particle results. The "disappearance" of the away-side 
correlation at high associated particle px, first observed 
for this kinematic range in the inclusive dihadron corre- 
lations in Ref. @, has a clear dependence on the trig- 
ger particle orientation. When the trigger particles move 
from (f> s ~ 0° to 90°, the path-length increases, so the 
away-side peak(s) become diminished. 

STAR has previously published dihadron correlations 
for in-plane (4> s < n/4) and out-of-plane (<j> s > n/4) trig- 
ger particles [22j . We sum up our correlation results from 
slices 1-3 and 4-6 to obtain the in-plane and out-of-plane 
correlations, respectively. We have also analyzed the data 
treating <p s < 7r/4 as a single slice to obtain the in-plane 
correlation and <f> s > 7r/4 for the out-of-plane correla- 
tion. The resultant correlation functions are consistent 
with those reported here that were summed from indi- 
vidual slices. Figure |S] compares results from this work 
to those in [22:]. The histograms show systematic uncer- 
tainties of the results from this work, while the shaded 
boxes show those of the results from 1221. The analysis 
reported here differs from that in Ref. [22j in two ways: 
(i) the average w 2 = («2{2,AS}+D2{4})/2 is used in this 
analysis while V2 ~ (^2{2} + «2{4})/2 containing signif- 
icant non-flow effect was used in Ref. [22j], (ii) the flow 
correlation is corrected up to 1)4 in this analysis while 
correction only up to t>2 was done in Ref. [22]. 

Figure [5] shows the in-plane and out-of-plane correla- 
tion results for two trigger pp ranges and two associated 

particle Py ranges for the 20-60% Au+Au collisions. 
The histograms show the systematic uncertainties due to 
flow subtraction; those due to zyam background normal- 
ization are shown as boxes in the legends. A difference 
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FIG. 8: (Color online) Comparison of dihadron correlation re- 
sults from this work (triangles) to those in Ref. [22| (circles) 
for (a) in-plane (<j> s < n/4) and (b) out-of-plane (<f) s > 
trigger particles. Data are from 20-60% Au+Au collisions. 
The trigger and associated particle pr ranges are 4 < pP < 
6 GeV/c and 2 < p ( p < 4 GeV/c respectively. Both the 
trigger and associated particles are restricted within \r}\ < 1. 
Error bars are statistical. Systematic uncertainties (including 
those due to anisotropic flow and due to background normal- 
ization from different zyam normalization ranges) are shown 
in histograms for results from this work and in shaded areas 
for results from Ref. [22[. The systematic uncertainty due 
to background deviations from ZYAM is not included for fair 
comparison. 



is observed between in-plane and out-of-plane dihadron 
correlations for both trigger Pj) ranges and both asso- 
ciated particle pp bins. The near-side correlated yield 
is larger for in-plane than out-of-plane triggers. As will 
be discussed in Sections [IV CI and IIV D[ the difference is 
due to the larger ridge contribution in-plane than out-of- 
plane, and the jet-like contributions are similar between 
in-plane and out-of-plane. A more significant difference 
is observed on the away side between in-plane and out- 
of-plane correlations. For in-plane trigger particles, the 
away-side correlations peak at A</> = 7r. For out-of-plane 
trigger particles, the away-side correlations are double- 
peaked. The double-peak structure is stronger for the 
lower trigger particle pP range. The away-side structure 
is studied in more detail in Sec. IIV Bl below. 



B. Discussion on the Away-Side Results 

In order to quantify the modification in the away-side 
correlation structure, we calculate the width of the dis- 
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FIG. 9: (Color online) Background-subtracted dihadron cor- 
relations with trigger particles in-plane {(f) s < tt/4) and out-of- 
plane (<f) s > 7r/4) in 20-60% Au+Au collisions. The results are 
for 3 < pP < 4 GeV/c (left panels) and 4 < pp < 6 GeV/c 
(right panels), and 1 < pip ' < 2 GeV/c (upper panels) and 
2 < pP < 3 GeV/c (lower panels). Both the trigger and 
associated particles are restricted within 77 j < 1. Error bars 
are statistical. Systematic uncertainties due to flow subtrac- 
tion are shown in histograms; those due to zyam background 
normalization are indicated by the vertical sizes of the filled 
and hollow boxes in the legends for in-plane and out-of-plane 
trigger particles, respectively. 
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We show in Fig. [TUTa) the RMS of the away-side corre- 
lation as a function of the trigger particle orientation 
(f> s for 20-60% Au+Au collisions. The associated parti- 
cle pr range is 1 < p^ < 2 GeV/c. Two trigger par- 
ticle pr ranges are shown: 3 < p>p < 4 GeV/c and 
4 < pip < 6 GeV/c. The RMS increases with increasing 
cj) s by nearly a factor of 2 from in-plane to out-of-plane. 
The distribution becomes more double-peaked as 4> s in- 
creases. No difference is observed between the two trigger 
pP selections. 

For comparison, the c?+Au results are indicated by the 
arrows in Fig. [TU^a) . As seen, the RMS in 20-60% Au+Au 
collisions from slices 1 and 2 is not much larger than in 
<i+Au. This may be consistent with the path-length ef- 
fect. However, we note that the correlation amplitudes in 
Au+Au collisions for the in-plane slices are larger than in 
d+Au collisions, as discussed below. This suggests that 
the away-side single peak in Au+Au and <i+Au collisions 
may come from different physics mechanisms. As will be 
discussed in Sec. IIV CI the near-side correlation for in- 
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of 4 < pj) < 6 GeV/c are qualitatively the same. 

The away-side double-peak structure observed in the 
inclusive dihadron correlation (i.e., without differenti- 
ating trigger particle azimuthal angles relative to the 
reaction plane) [To| has stimulated much interest [l7l - 
[2pj . The recent three- particle jet-like correlation studies 
have provided evidence of conical emission of hadrons 
correlated with the high pr trigger particles [21] . To 
study conical emission in more detail, we show in Fig. [TT] 
the average correlation amplitude on the away side in 
the 7r-region (\A<j> — tt\ < 0.39) and in the cone-region 
(0.81 < \A<f> - tt| < 1.59) as a function of (j) s in 20-60% 
Au+Au collisions. Two trigger pp) selections are shown; 
no significant difference is observed. The amplitudes in 
the tt- region drop with increasing cj> s , from a value larger 
than that in d+Au (as indicated by the arrows) to a 
value significantly smaller than that in d+Au. On the 
other hand, the cone-region amplitude seems rather con- 
stant with 4> s in Au+Au collisions. The amplitude in 
cone-region is significantly stronger than that in d+Au 
collisions for both trigger particle pp selections. This 
suggests that conical emission may also be present for 
in-plane trigger particles. See discussion in Sec. IIVDI 



FIG. 10: (Color online) (a) The away-side RMS of the dihadron 
correlation function versus the trigger particle azimuth rela- 
tive to the event plane, 4> s = \<j>t — 4>ep\, in 20-60% Au+Au 



collisions for 1 < p^ 1 < 2 GeV/c. Two trigger pP selec- 
tions are shown: 3 < pP < 4 GeV/c (solid circles) and 
4 < pP < 6 GeV/c (hollow squares), (b) The away-side RMS 
for slice 1 (hollow triangles) and slice 6 (solid triangles) versus 
the associated particle p)p ' in 20-60% Au+Au collisions. The 

trigger particle pr range is 3 < pp < 4 GeV/c. Error bars 
are statistical. The curves indicate systematic uncertainties 
due to flow subtraction, solid curves for the solid data points 
and dashed curves for the hollow data points. The system- 
atic uncertainty due to zyam background normalization is not 
shown. The corresponding d+Au results are indicated by the 
arrows (solid arrow for 3 < Pj? < 4 GeV/c and hollow arrow 
for 4 < pP < 6 GeV/c) in the upper panel and by the lower 
solid line connecting error bars in the lower panel. 



plane trigger particles has a large contribution from the 
ridge, and it is likely that there is an accompanying back- 
to-back ridge on the away side. 

Figure ITTjT b) shows the RMS as a function of the asso- 
ciated particle pp* for slices 1 and 6 in 20-60% centrality. 
The RMS remains constant for slice 1, and is not much 
broader than the d+Au result for all measured pp bins. 

The RMS for slice 6 increases with pp\ The double- 
peak structure is the strongest when the trigger particle 
is perpendicular to the reaction plane and the associated 
particle is hardest. Results for other slices vary smoothly 
between slices 1 and 6. The features for trigger particles 
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FIG. 11: (Color online) The away-side dihadron correlation 
amplitudes in the 7r-region (\A<f> — 7r| < 0.39) and the cone- 
region (0.81 < \A<f> — 7r| < 1.59) as a function of the trigger 
particle azimuth relative to the event plane, cf> s = \4>t— V'Epj, in 
20-60% Au+Au collisions. Statistical errors are smaller than 
the symbol size. The curves indicate systematic uncertainties 
due to flow subtraction, and the brackets indicate those due to 
zyam background normalization. Both trigger pr ranges are 
shown: (a) 3 < pP < 4 GeV/c and (b) 4 < p { p < 6 GeV/c. 
The associated particle pr range is 1 < p^? < 2 GeV/c. The 
d+Au results in the tt- and cone-regions are indicated by the 
solid and hollow arrows, respectively. 
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Comparison of the relative amplitudes in the 7r-region 
and the conical emission region shown in Fig. [TT] again 
reveals the degree of the double-peak structure. In order 
to study the px dependence of the relative amplitudes, 
Fig. HHa) and (b) show the amplitude ratios of 7r-region 
to cone- region in slices 1 and 6, respectively. The ampli- 
tude ratio in slice 1 increases with p^ ' for the higher pj) 
trigger particles. The trend is not much different from 
that observed in d+Au collisions (shown by the black 
line). The increasing trend suggests that for in-plane 
trigger particles the away-side correlation is dominated 
by physics mechanisms other than conical emission, such 
as punch-though jets and/or back-to-back ridge. The in- 
creasing trend may also be present for the lower p^) trig- 
gers, but the systematic uncertainty in this analysis pre- 
vents a firm conclusion. On the other hand, for slice 6 the 

amplitude ratio decreases with p^ . The away-side jet- 
like correlation at Acf) = ir is essentially diminished; what 
remains are conical emission hadrons. It is also worth to 
note that the away-side amplitude ratio for out-of-plane 
trigger particles (lower panel of Fig. [TT]) is significantly 
smaller than for in-plane trigger particles (upper panel of 
Fig. HT|) . This is again the consequence of the significant 
away-side broadening from in-plane to out-of-plane. 



C. Discussion on the Near-Side Results 

Recall that in Fig. [7] we observe a significant change 
in the near-side peak amplitude. The near-side ampli- 
tude drops with increasing <j) s . For the 20-60% centrality, 
the amplitude at large (f> s is not much different from the 
d+Au result, perhaps indicating minimal medium modifi- 
cation. On the other hand, the amplitude at small <p s ap- 
pears significantly larger than in d+Au suggesting large 
medium effect. This might be counterintuitive at the first 
look. Due to jet- quenching, the near-side jets predomi- 
nately emerge outward from the surface of the medium, 
so variation in the medium thickness the near-side jets 
traverse between in-plane and out-of-plane directions is 
not naively expected. 

It has been shown by the inclusive dihadron correla- 
tion that the near-side correlation strength is enhanced in 
Au+Au with respect to p+p and d+Au collisions (Tolll^. 
and the enhancement is mainly due to the large contri- 
bution from the ridge [10, fll| . In order to investigate 
the underlying physics mechanism for the near-side struc- 
ture change with trigger particle orientation, we separate 
contributions from the ridge and the jet-like component 
by analyzing the correlation data in two different A?y 
regions [ll|: IA77I > 0.7 where the ridge is the domi- 
nant contributor and |A?y| < 0.7 where both the ridge 
and jet-like correlations contribute. Figure Q2] (upper 
panel) shows the background-subtracted dihadron cor- 
relation function from \Ar)\ > 0.7 for trigger and asso- 
ciated particle px ranges of 3 < Pj) < 4 GeV/c and 
1 < p^ ] < 2 GeV/c in 20-60% Au+Au collisions. (The 
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FIG. 12: (Color online) Associated particle dependence 
of the ratio of away-side dihadron correlation amplitude in 
the 7r-region (|A0 — tv\ < 0.39) to that in the cone-region 
(0.81 < |A0 - 7r| < 1.59). Two 4> s = \4> t - ip E p\ slices are 
shown: (a) 0° < 4> s < 15° and (b) 75° < 4> s < 90°. The data 
are from minimum-bias 20-60% Au+Au collisions. Both trig- 
ger pt ranges of 3 < p^! < 4 GeV/c and 4 < p)p < 6 GeV/c 
are shown. Error bars are statistical. The curves indicate 
systematic uncertainties due to flow subtraction. The sys- 
tematic uncertainty due to zyam background normalization 
is not shown. The d+Au results are indicated by the open 
triangles connected by the line in (a), where the error bars 
are statistical. 



|Afy| > 0.7 correlation functions for all kinematic ranges 
are presented in Appendix [A] in Figs. [55] and |2"9"]) . The 
near-side correlation for | A77I > 0.7 is due to the ridge be- 
cause the jet-like contribution is mostly confined within 
IA77I < 0.7. The ridge correlation shows a significant 
drop with increasing tf> s . The ridge contribution is close 
to zero for trigger particles perpendicular to the reaction 
plane in the 20-60% centrality bin. 

It has been recently suggested by the NeXSPheRIO 
model [64T466T ] that initial energy density fluctuations (hot 
spots) with subsequent hydro evolution may generate a 
near-side ridge and a double-peak correlation on the away 
side. The physics mechanism appears to be side-splashes 
of particles by the hot spot on the surface resulting in 
two peaks in the single particle azimuthal distribution 
event-by-event separated by about two radians [HI, H|| . 
These two-peaked single particle distributions produce 
two-particle correlations of a near-side ridge and an away- 
side double-peak. The near-side ridge and the away-side 
double-peak are due to the same physics, and the near- 
side ridge amplitude should be larger (by a factor of two) 
than each of the two away-side peaks. This relative am- 
plitude is a unique feature of the NeXSPheRIO model 



20 



0°-15° 



15°-30° 



30°-45° 



45°-60° 



60°-75° 



75°-90° 



-e- 0.3 

< 0.2 

Z 01 

T3 o 




-0.1 



20-60% Au+Au, 3<pf <4 GeV/c, 1<pf ><2 GeV/c 



i ■ ■■ i ■ ■■ i 



i ■ ■■ i ■ ■■ i 



i ■ ■■ i ■ 1 1 1 



i ■ 1 1 1 ■ ■■ i 



i 1 ■■ i ■ ■■ i 



a 

ID 
CD 



024 024 024 024 024 024 

A<|) = <j) - <|) [rad] 

FIG. 13: (Color online) Upper panels: background-subtracted dihadron correlations with trigger particles in six slices of 
azimuthal angle relative to the event plane, cj> s — \4>t — V'epI, with a cut on the trigger-associated pseudo-rapidity difference 
of \Ari\ > 0.7. Statistical errors are smaller than the symbol size. Systematic uncertainties are shown in the black histograms 
due to flow subtraction and in the horizontal shaded band around zero due to ZYAM background normalization. The near-side 
correlation is due to the ridge. Lower panels: The difference between raw dihadron correlations in \Ar)\ < 0.7 and \Arj\ > 0.7, 
after multiplying a coefficient onto the later such that the resultant difference is zero on average on the away side in the range 
|A</> — ir\ < 1. This correlation represents the jet-like component of the dihadron correlations. Statistical errors are smaller than 
the symbol size. Systematic uncertainties are small. The results are for 3 < < 4 GeV/c and 1 < < 2 GeV/c in 20-60% 
Au+Au collisions. In both panels the corresponding inclusive dihadron correlations from d+Au collisions (thick histograms) 
are superimposed for comparison. 



because of the topology of particle distributions in the 
model [HI, [66| . This feature is not observed in the out- 
of-plane large Ar) correlations in data. 

It has also been recently suggested based on the AMPT 
and UrQMD models [47ll49j that there may be large tri- 
angularity in the initial collision geometry event-by-event 
and those initial geometry fluctuations could produce a 
triangular anisotropy (triangular flow) in the final mo- 
mentum space. Such triangular flow would result in three 
peaks at A0 = 0, 2ir/3, and 4ir/3 in the two-particle cor- 
relation, which appear to be qualitively consistent with 
the inclusive dihadron correlation data integrated over 
all reaction plane directions [l3|, [TJ]. The authors of 
Ref. [13, HH thus attribute the ridge and the away-side 
double-peak to possibly triangular flow. Because the mi- 
nor axis direction of the initial fluctuating triangular ge- 
ometry is random with respect to the reaction plane or 
the participant plane, the three-peak structure in two- 
particle correlation from triangular flow should be inde- 
pendent of <p s . However, the near-side peak of our di- 
hadron correlation data decreases with increasing <p a and 
is consistent with zero at large Ar] with trigger particles 
out-of-plane, as shown in the upper panel of Fig. [UU as 
well as in Figs. [28l and [29l These results suggest that tri- 
angular flow is not a major contribution to the measured 
dihadron correlations with high px trigger particles. 

We may estimate the effect of triangular flow fluctua- 
tions. In our 20-60% centrality range, model studies (47i — 
l49j as well as inclusive two-particle correlation measure- 
ments [50| indicate v\jv\ ~ 0.1 for our trigger and asso- 
ciated pt ranges. Using the V2 parameters in Table U and 
the background magnitudes in Table IIII1 we obtain the 



i>3 effect on the dihadron correlations at |Ar/| > 0.7 with 
p ( V =3-4 GeV/c to be ±2Bv ( 3 a) vf « ±0.023, ±0.028, 

±0.011, ±0.004, ±0.002, for p^ ] =0.15-0.5, 0.5-1, 1-1.5, 
1.5-2, and 2-3 GeV/c, respectively. Compared to the cor- 
relation strengths in Fig. [2H1 these account for only about 
10% of the measured correlations for the respective p^ 
bins. 

In order to make a quantitative estimate of the 
t>3 effect on our dihadron correlations, we include a 
±2v^v^ cos(3A0) term in Eq. ([1]) for our flow back- 
ground. We then apply the ZYAM normalization proce- 
dure for background subtraction. Figure Q3] shows the 
obtained dihadron correlation functions for \Ar)\ > 0.7 
and 3 < pJ^ < 4 GeV/c. We used the estimated v\jv\ of 
0.1 (upper panels), as well as values twice (middle panels) 
and three times larger (lower panels). One notices that 
the near-side correlations for out-of-plane trigger parti- 
cles become concaved with subtraction of v\jv\ = 0.2 
and 0.3, indicating an oversubtraction of the triangular 
flow effect. As seen from Fig. [131 the qualitative features 
of the correlation functions are unchanged from those in 
the upper panels of Fig. [T3J The away-side double-peak 
structure out-of-plane remains prominent, even with the 
significantly larger i>3 subtraction. The decreasing trend 
of the ridge magnitude from in-plane to out-of-plane is 
unaffected because a constant contribution over <f) s is 
subtracted. Our results discussed later in the paper are 
also not affected qualitatively. These demonstrate that 
the main features of the measured near-side ridge and 
away-side double-peak in the dihadron correlations with 
high pt trigger particles, whether or not integrated over 
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FIG. 14: (Color online) Effect of triangular flow fluctuations on dihadron correlations at |A?7| > 0.7 for 3 < p§ < 4 GeV/c and 
1 < ' < 2 GeV/c in 20-60% Au+Au collisions. The effect is independent of the trigger particle azimuth relative to the event 
plane, cj> s . The flow systematic uncertainties are shown in the black histograms. Three values of triangular flow are shown; 
«!/i>2 = 0.1 in the upper panels (representing the best estimate of 1)3), 0.2 in the middle panels, and 0.3 in the lower panels. 
The solid curves represent the subtracted triangular flow effect. 



<f> s , are not due to the possible triangular flow fluctua- 
tions, but other physics mechanisms. 

The near-side ridge correlation at large A77, after two- 
particle Arj acceptance correction, was found to be nearly 
uniform in A77 [111 ]. If the ridge is uniform over the en- 



tire measured A77 range, then the ridge can be readily 
subtracted by taking the difference between the raw (not 
background-subtracted) correlations from the small and 
large A77 regions as 
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cWjet-li 
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dA^ dA V 



(12) 



The coefficient A accounts for the A77 acceptance differ- 
ence between |A?y| < 0.7 and \Ar]\ > 0.7, and can be 
easily obtained from the acceptance ratio of the two Arj 
regions. It can also be obtained by requiring the away 
side of the resultant average correlation magnitude to 



be zero because the away-side correlation (after Arj ac- 
ceptance correction) is also uniform within the measured 
A77 range in the TPC [Tij. We use the latter method 
to obtain A such that the resultant away-side average 
correlation signal within \A<j> — n\ < 1 is zero, namely 
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(13) 



The obtained coefficient is approximately A ~ 1.45. The 
resultant difference by Eq. (fj"2"j) represents the dihadron 
correlation of the near-side jet-like component under the 



assumption that the near-side ridge is uniform in Arj 
within our measured range. The Acj) correlation of the 
jet-like component obtained by Eq. (fT2"]) is free of large 
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systematic uncertainties because the anisotropic flow is 
cancelled in the difference. 

The obtained Acf> correlation of the jet-like component 
is shown in the lower panel of Fig. 1131 The corresponding 
d+Au result is superimposed on the figure. The A<j> cor- 
relation of the jet-like component is approximately inde- 
pendent of the trigger particle orientation, in contrast to 
the ridge component shown in the upper panel of Fig. 1131 
The near-side jet-like correlations are consistent between 

I 



Ridge yield = 




We extract the jet-like yield in |A»j| < 0.7 and \A<j>\ < 1 
from the correlations of the jet-like component, such as 
those in the lower panel of Fig. [T3l by 

1 f 1 f - 7 d 2 N 

Jet-like yield = — - / dAcj) / — — dA v . 

iVtrig J_ 1 J_ 7 dAcpdAr] 

(15) 

Note the Ar] acceptance is not corrected in the A<j> cor- 
relations of the ridge or the jet-like component; hence, 
neither are the extracted corresponding yields. The ex- 
tracted ridge and jet-like yields are shown in Fig. [15] as 
a function of <f> s in the 20-60% centrality bin. The boxes 
indicate the total systematic uncertainty; the individual 
sources of systematic uncertainties and their correlations 
have been discussed earlier in Sec. IIII1 As seen from 
Fig. 1151 the jet-like yield is approximately independent 
of 4> s in Au+Au collisions, and consistent with the d+Au 
data. The ridge yield in Au+Au collisions at small cf> s 
(in-plane) is significant, but it decreases quickly with in- 
creasing 4> s . The ridge yield at large 4> s (out-of-plane) is 
consistent with zero. The ridge is dominated by events 
where trigger particles are within approximately 45° of 
the event plane. 

The trend of decreasing ridge amplitude with increas- 
ing <j) a is seen in all measured p^ bins. To quantify 
this, we show in Fig. ITBTa) and (b) the p^ depen- 
dence of the ratio of ridge yield in 75° < <j> s < 90° and 
30° < (j) s < 45°, respectively, to that in 0° < cj) s < 15°. 
Both trigger particle Pj) selections are shown. The sys- 
tematic uncertainties, shown for 3 < p^) < 4 GeV/c in 
the boxes, have taken into account correlations among 
the different sources of systematic uncertainties. Within 
the systematic uncertainty there is no observable differ- 
ence between the two Pj) selections. The ridge ratios 
from different <f> s slices appear to be independent of p^ . 

The ridge decreases with 4> s universally for all pi? . The 
ridge yield out-of-plane is consistent with zero at all as- 
sociated particle p>p' for both the trigger particle pip 
selections. 



d+Au and Au+Au collisions. The A(j> correlation func- 
tions of the jet-like component for all trigger and associ- 
ated particle px ranges are presented in Appendix [X] in 
Figs. GUI and EU 

To quantify the near-side modification, we study the 
ridge and jet-like yields as a function of 4> s . We extract 
the ridge yield in \Arj\ > 0.7 and \A(/)\ < 1 from the zyam 
background-subtracted correlations, such as those in the 
upper panel of Fig. 1131 by 



dA^, dAT1 + L dAHA^ dA V ■ (14) 
I 

Motivated by the preliminary version of our data, Chiu 
and Hwa [(57j suggested that alignment between jet prop- 
agation and medium flow direction, likely to be found 
for in-plane trigger particles, may be responsible for the 
ridge; radiated gluons (within a small angle of the par- 
ton direction) become thermalizcd with the medium and 
combine with medium partons to form the ridge when 
they are aligned in the same direction. We note that this 
correlated emission of ridge particles with the medium 
flow direction may be rather general, not necessarily re- 
stricted to the recombination of radiated and medium 
gluons. For instance, it is possible that initial fluctua- 
tions of color flux tubes together with the stronger in- 
plane transverse flow can produce similar effects [64|, |68| - 
|72I ]. We discuss this color flux tube fluctuation model 
further in Sec. llVDl 

There is strong experimental evidence suggesting that 
the jet-like component and the r idg e are produced by 
different physics mechanisms [HI, |4l[ , thus their px de- 
pendences are expected to be different. To quantitatively 
study this, we show in Fig. RTOtc) the ratio of the ridge 
yield to the jet-like yield for 0° < cf) s < 15°. Again, 
the systematic uncertainties shown in boxes have already 
taken into account correlations among different sources 
of systematics. Within the systematic uncertainties, the 
ridge over jet-like component ratio appears to be con- 
stant over the measured j?y . This may suggest, on the 
contrary to the other findings, that the ridge and the jet- 
like component may be of the same origin. However, it is 
possible that differences in the p^ spectra of the jet-like 
and the ridge component are small for our trigger p^ 
ranges compared to our systematic uncertainties. The 
Pt spectra of the jet-like component and the ridge will 
be further discussed below. 



D. Connections between Near- and Away-Side 

We have observed that the away-side amplitude in 
the 7r-region decreases strongly with increasing <j) Sl so 
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(a) 3<p™<4 GeV/c, 1<p^<2 GeV/c 

• jet(A(|>|<1.0, |AT]|<0.7) 
■ ridge (A<H<1.0, |Ati|>0.7) 
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^ S =l^ t -¥ EP I [deg] 

FIG. 15: (Color online) The near-side jet-like and ridge yields 
as a function of the trigger particle azimuth relative to the 
event plane, <j> s = \<f>t — ipEp\- The results are from 20-60% 
Au+Au collisions. Two trigger ranges are shown: (a) 



3 < p%> < 4 GeV/c and (b) 4 < p^> < 6 GeV/c. The associ- 
ated particle p^ range is 1 < p^f ' < 2 GeV/c. The jet-like 
yield is from \Acj>\ < 1 and \Ar)\ < 0.7 and the ridge yield 
is from \Acj>\ < 1 and |A?7 > 0.7. Error bars are statistical. 
The systematic uncertainties are shown by the boxes. For 
the ridge yield they include those from anisotropic flow (in- 
dicated by the curves) and zyam background normalization. 
The systematic uncertainties on the jet-like component are 
due to leakage of jet-like correlations out to |A^| > 0.7 and 
the assumption that the ridge is uniform in Arj. The d+Au 
results in the jet and ridge regions are indicated by the filled 
and hallow arrows, respectively. 
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FIG. 16: (Color online) (a) Ratio of the ridge yield from 75° < 
<t>s < 90° to that from 0° < 4> s < 15°. (b) Ratio of the ridge 
yield from 30° < 4> a < 45° to that from 0° < <t> a < 15°. (c) 
Ratio of the ridge yield to the jet-like yield from 0° < 4>s < 
15°. The ridge yield is from |A^| < 1 and \Ar/\ > 0.7 and 
the jet-like yield is from |A<^| < 1 and \Arj\ < 0.7. Data are 
from 20-60% Au+Au collisions. Both trigger ranges of 
3 < p$ < 4 GeV/c and 4 < p^} < 6 GeV/c are shown. Error 
bars are statistical. Boxes indicate systematic uncertainties 
on the 3 < pP < 4 GeV/c data. 



does the near-side ridge amplitude. We have also ob- 
served that both the away-side amplitude in the cone- 
region and the near-side jet-like amplitude remain ap- 
proximately constant with cj> s . This raises the question 
whether the near-side and the away-side are connected, 
or stem from the same physics origin, even though high 
Pt trigger particles strongly bias the near-side towards 
surface emission. In order to gain further insights, we 
study the near- and away-side correlation properties to- 
gether as a function of <j) s and . 

Figure [T7l shows the average correlation amplitudes of 
the away-side 7r-region and cone-region and the near-side 
ridge and jet-like component. The averages are taken 
within the same window size of ±0.39. The ridge ampli- 
tude is scaled by a factor of 2.45, which is approximately 
the acceptance factor to scale |Ar;| > 0.7 to the entire 
Arj range assuming a uniform ridge. The jet-like ampli- 
tude and the cone-region amplitude have a similar depen- 
dence on <j) s . The similarity suggests that the near-side 
jet-like component and the away-side cone-region might 



be closely related. This is expected when the away-side 
parton loses a significant fraction of its energy into the 
conical emission of particles [201 ] . 

Figure [17] also shows that the ridge amplitude and the 
away-side 7r-region amplitude have a similar dependence 
on <j) s . The magnitudes are also similar between the ridge 
and the 7r-region. This is especially true for the lower pj) 
range. On the other hand, the jet-like and cone-region 
amplitudes have a rather different dependence on <fi a than 
the ridge and 7r-region amplitudes. This suggests that 
the near-side ridge and the away-side 7r-region may be 
connected. On the contrary, they seem not connected to 
the jet-like component or to the component in the cone- 
region. 

There are many other experimental evidences suggest- 
ing that the ridge and the jet-like component may be 
unrelated in physics despite of the apparent correlation 
between the ridge and the high px trigger particle. For 
example, three-particle correlations suggest that the pro- 
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TABLE IV: Inverse slope parameter T (MeV/c) from an exponential fit to the associated particle pr spectra of correlated 
amplitudes in different A<f> regions: dN / [prdpr) oc exp(— pr/T). Systematic uncertainties for the jet-like spectra and the 
(•7r-region — ridge) spectra are small. Statistical and systematic uncertainties on inclusive charged hadron spectrum are both 



negligible. 




0° < 4 
3 < pf < 4 GeV/c 


> s < 45° 

4 < < 6 GeV/c 


45° < ( 
3 < pip < 4 GeV/c 


b s < 90° 
4 < p ( p < 6 GeV/c 


Jet-like 
Cone-region 
7r-region 
Ridge 
7r-region — Ridge 


465 ± 6(stat) 
343 ±2(stat)+g 4 (syst) 
357 ± l(stat)+* (syst) 
432 ±3(stat) + *| (syst) 


518 ± 13(stat) 
316±4(stat)+g 1 (syst) 
355 ±3(stat)^ 2 (syst) 
417±8(stat)+^(syst) 
239 ± 14(stat) 


460 ± 6(stat) 

341 ± 2(stat)+f (syst) 

238±3(stat)+^ (syst) 


518 ± 14(stat) 

333 ±4(stat)+7 2 (syst) 

288 ±7(stat)+^ (syst) 

252 ± 14(stat) 


Inclusive charged hadron 
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duction of the jet-like component and the production of 
the ridge are uncorrelated [l2[ . The particle composition 
of the ridge was found to be similar to that of the bulk 
medium 4l|. The ridge magnitude was observed to be 
rather independent of the trigger particle pJp, persist- 
ing to very large Pj) [HI where jets are almost the sole 
source of those large Pj) trigger particles. The parent 

parton energies triggered by the wide range pp trigger 
particles vary greatly, and yet the ridge is independent of 
pip . This, again, suggests that the ridge and the jet-like 
component may be unrelated. 

It has been recently suggested that the ridge may be 
generated by fluctuations of color flux tubes stretched 
between the colliding nuclei at the initial time of con- 
tact [63, [ZQHzli • The ridge particles from the color flux 
tubes near the surface of the collision zone are boosted 
radially by the medium expansion, becoming correlated 
in relative azimuth. If the ridge is indeed due to color flux 
tube fluctuations, i.e. entirely from the medium without 
connection to high px trigger particles, then the meaning 
of "near side" as defined by the high pr trigger particle 
bears no significance to the ridge. In such a case, there 
ought to exist a ridge partner on the away side due to 
symmetry, i.e. a back-to-back ridge. In addition, it is 
conceivable that the ridge would be stronger along the re- 
action plane direction where both the flux tube strength 
and the medium flow are stronger. This would natuar- 
ally explain our observation that the ridge decreases from 
in-plane to out-of-plane and the ridge amplitude and the 
7r-region amplitude trace other other. The trigger parti- 
cles in the in-plane direction happen to have ridge par- 
ticles associated within a narrow A<f> region (near-side), 
while those trigger particles out-of-plane cannot acciden- 
tally pick up ridge particles to be within a narrow A<f> 
azimuth. In fact, the above mechanism where the ridge 
particles are aligned with the trigger particle in azimuth 
is similar to the ridge formation mechanism proposed in 
the Correlated Emission Model [g^]; however, the under- 
lying physics is quite different. 

Examining the px dependences of the different cor- 
relation components can give further insights into the 



physics mechanisms responsible for their formation. We 
show in Fig. [TS] the pp ' spectra of the average correla- 



tion amplitudes from various A<fi regions, 



dN 
PrdpT ' 



The 



upper panels show results for in-plane trigger particles, 
0° < 4> s < 45°. Four Atfi regions are shown: the ir- region, 
the cone-region, the jet, and the ridge. The lower panels 
show results for out-of-plane particles, 45° < <f) s < 90°. 
The ridge, which is consistent with zero, is not shown 
in the lower panel for clarity. Both pip selections are 
shown, 3 < pj) < 4 GeV/c in the left panels and 

4 < pP < 5 GeV/c in the right panels. Note that the 
upper panels and the lower panels have the same order 
of magnitude span in their coordinates, so the spectral 
shapes can be readily compared. 

To investigate the spectral shapes quantitatively, we 
fit an exponential function to the spectra. The inverse 
slopes from the fits are tabulated in Table IIVI The sys- 
tematic uncertainties of the fitted inverse slope param- 
eters have already taken into account the correlations 
in the various sources of systematic uncertainties of the 
spectra. As expected, the jet-like spectra are harder for 
the higher trigger pp) range. The difference in the jet-like 
spectra between in-plane and out-of-plane is small. On 
the other hand, the cone spectra do not seem to depend 
on trigger particle pp , nor trigger particle orientation 
relative to the reaction plane. The cone-region appears 
to be universal. In addition, the cone spectra are signif- 
icantly softer than the jet-like spectra, suggesting a dif- 
ferent production mechanisms for the near-side jct-likc 
correlated hadrons and the away-side conical emission 
hadrons. Yet, the hadron yields in the jet-like correla- 
tion region and in the cone-region appear to trace each 
other. This would be a natural consequence if the away- 
side parton, in rough energy balance with the near-side 
jet, loses most of its energy to conical flow [jq . 

It is interesting to note that the 7r-region hadrons 
are similar to the conical emission hadrons in the in- 
plane direction; however, in the out-of-plane region they 
are softer. In fact, the out-of-plane hadrons in the tt- 
region are not much different from the inclusive hadrons 
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FIG. 17: (Color online) Average correlation amplitude as a 
function of (j> a = \4>t — V'epI for the near-side jet-like compo- 
nent (\A<j}\ < 1, |Ar?| < 0.7), the cone-region (|A0 — 7r±1.2| < 
0.39, \r)\ < 1), the 7r-region (|A0 - vr| < 0.39, \<r)\ < 1), and 
the acceptance-scaled near-side ridge (|A0| < 1, \Arj\ > 0.7). 
Data are from 20-60% Au+Au collisions. Both trigger pP 
selections are shown: (a) 3 < pip < 4 GeV/c and (b) 
4 < pP < 6 GeV/c. The associated particle pr range is 
1 < pi^ < 2 GeV/c for both panels. Error bars are statistical. 
Boxes indicate systematic uncertainties due to anisotropic 
flow. The systematic uncertainties due to zyam background 
normalization, common to the 7r-region, cone region and ridge 
amplitudes, are not shown. 



in their pt distributions. In the scenario of only jet- 
quenching, the in-plane away-side partons do not have 
enough medium to interact in 20-60% Au+Au collisions 
to completely wash out their identity. On the other hand, 
the out-of-plane away-side partons have a longer path- 
length and the lost energy appears to have equilibrated 
with the medium, a result found in the inclusive dihadron 
correlation in central collisions [Toj . 

Surprisingly, the ridge particles are relatively hard, 
not much softer than the jet-like particle spectra (see 
Fig. HH and Table. [IVj) . Yet, the <p s dependence of the 
ridge yield is completely different from that of the jet- 
like yield. We note that the ridge spectrum measured at 
larger pip > 2 GeV/c is significantly softer than the jet- 
like hadron spectrum also at large pip suggesting 
that the ridge might be related to the medium. If the 
ridge comes from the medium, then our result implies 
that it is not a simple uniform share of medium parti- 

» 



particles are harder than the bulk medium particles. 

For the associated particle pr range of 1 < pi^ ' < 
2 GeV/c shown in Fig. [TTl the 7r-region amplitude is 

slightly smaller than the ridge amplitude for 3 < pip < 

4 GeV/c. For lower Pj?\ the 7r-region amplitude is gen- 
erally larger than the ridge amplitude. For the higher 
trigger px range of 4 < pip < 6 GeV/c, there appears an 
excess of particles in the away-side 7r-region over those in 
the near-side ridge for all pip bins. The excess appears to 
be insensitive to <p s . Experimentally, it is interesting to 
examine the px dependence of this excess by taking the 
difference between the away-side 7r-region and the near- 
side ridge. This difference is rather robust because all the 
systematic uncertainties cancel. The difference (excess 
of particles in the away-side 7r-region over the near-side 
ridge) is shown as diamonds in Fig. [THJ It is remarkable 
to note that the pr spectra of those excess particles are 
rather similar between in-plane and out-of-plane in terms 
of their inverse slopes. Direct comparison is made in the 
lower panels where the excess particle spectra in the up- 
per panels are superimposed in the corresponding lower 
panel as dotted lines. The agreement is excellent for the 
trigger particle 4 < pip < 6 GeV/c range while the statis- 
tics are insufficient at large pP for the lower trigger pip 
range. We found the agreement persists in each individ- 
ual 4> s slice. Those excess particles have rather soft Pt's, 
similar to the inclusive charged hadrons. This is already 
evident from the ridge and 7r-region spectra; the away- 
side 7r-region spectra are softer than the ridge spectra. 
If the ridge is generated by fluctuating color flux tubes 
and is back-to-back [H, [7(J|72[ , then the excess particles 
in the away-side 7r-region must come from other physics 
mechanisms. One such mechanism is punch-through jets. 
However, it is counterintuitive to have a much softer spec- 
trum for punch-through jet-like particles, as well as the 
agreement between in-plane and out-of-plane directions. 
Another mechanism is statistical global momentum con- 
servation to balance the extra momentum carried by the 
ridge particles (because they are harder than the particles 
in the 7r-region). However, one may expect a somewhat 
harder spectrum for the recoil from statistical global mo- 
mentum conservation than the inclusive spectrum [73J . 



E. Properties of the Correlation Peaks 

To characterize the structure of the correlation func- 
tions, we fit the large-A?7 azimuthal correlations with 
two away-side Gaussian peaks symmetric about Acf> = tt 
and two ridge Gaussians (at A(f> = and ir) with iden- 
tical widths. We allow the ridge Gaussian magnitudes 
to vary independently because physics mechanisms other 
than the back-to-back ridge can also contribute to the 
7r-region as discussed earlier. The fit results are shown 
by the curves in Fig. [15] for 3 < pip < 4 GeV/c and 

(a) 



cles at our measured relatively low p T because the ridge 1 < p T < 2 GeV/c as an example. The away-side to 
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FIG. 18: (Color online) Correlated particle pr spectra in different A<f> regions: the near-side jet-like component (\A<f>\ < 1, 
A?;| < 0.7), the cone-region (|A</> - 7T ± 1.2| < 0.39, \rj\ < 1), the 7r-region (|A(/> - 7r| < 0.39, |r?| < 1), and the near-side 
ridge (|A<^> < 1, \Arj\ > 0.7). The ridge amplitude is scaled by the two-particle A77 acceptance ratio of approximately 2.45. 
The data are from minimum-bias 20-60% Au+Au collisions. Error bars are statistical. Systematic uncertainties (including 
those on anisotropic flow and on zyam background normalization) are shown as boxes for the cone-region, the 7r-region, and 
the ridge spectra. Two trigger pip ranges and two (f> s regions are shown: (a) 3 < pp 1 < 4 GeV/c and 0° < (fa < 45°, (b) 
4 < pP < 6 GeV/c and 0° < <j> a < 45°, (c) 3 < p£ 5 < 4 GeV/c and 45° < <j> s < 90°, and (d) 4 < p$> < 6 GeV/c and 
45° < 4> s < 90°. The dot-dashed curve is the inclusive charged hadron spectrum with an arbitrary normalization. All other 
lines are to guide the eye. The dotted curves in (d) is replicate of the dashed curve in (b). 



near-side ridge ratio for 3 < pj, < 4 GeV/c is generally 

larger than unity at low p^ but becomes smaller than 

unity at large \ For 4 < p^p < 6 GeV/c the fit error 
is too large to draw a firm conclusion. 

We study the peak positions and the Gaussian widths 
of the various components in the dihadron correlation ob- 
tained from our four-Gaussian fit. The Gaussian widths 
are shown in Fig. l2"0Ta) as a function of 4> s for the asso- 
ciated particle pft range of 1 < p^ < 2 GeV/c. Also 
shown are the jet-like peak widths in A77 fit to the near- 
side At? correlation functions (\A<fi\ < 1), as in Fig. [5] 
The jet-like peak widths as a function of <j) s are flat and 
are consistent between A(f> and A77 for this associated 
1 < pfi < 2 GeV/c bin. This again indicates that the 
jet-like correlation component is independent of the ori- 
entation of the trigger particle. The width for the con- 
ical emission peaks increases somewhat with increasing 
<j) s . This may be due to a more significant influence of 
medium flow on the conical emission angle out-of-plane 
than in-plane as discussed below. The ridge width de- 



creases somewhat with increasing <f> s , while the magni- 
tude of the ridge decreases significantly with <j) s , as shown 
earlier. 

Figure [2Uf b) shows the peak Gaussian widths as a func- 
tion of p^ 1 for integrated 4> s . The jet- like width in A</> 

decreases with increasing associated particle p^ , consis- 
tent with expectations for jet fragmentation. The A77 
width of the jet-like component generally agrees with the 
A<j) width for p^ > 1 GeV/c, but appears narrower than 
the A<f> width for lower p^ . The conical emission peak 

width does not vary significantly with p^r , and is wider 
than the near-side jet-like peak. It is interesting to note 
that the ridge peak width does not depend on associ- 
ated particle p^ . This is in contrast to the A<f> width of 
the jet-like peak. This may be taken as a confirmation 
that the ridge and the jet-like component, both on the 
near side of the trigger particle, may come from rather 
different physics mechanisms. 

For comparison, we fit the near-side jet-like A<f> corre- 
lation in minimum bias <i+Au collisions [which was also 
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FIG. 19: (Color online) Fit to the A</> correlation functions at |A?j| > 0.7 in six slices of <f) a = \4>t — V'epI by four Gaussians 
(curves): a near-side Gaussian at A<f> = corresponding to the ridge, an away-side Gaussian at A<f> = 7r with identical width to 
the near-side Gaussian but with varying amplitude, and two identical away-side Gaussians symmetric about A(j> — it. Data are 
from 20-60% Au+Au collisions, as same as those in the upper panel of Fig. [13] The trigger and associated particle pr ranges 
are 3 < pfj? < 4 GeV/c and 1 < Pj^ < 2 GeV/c, respectively. The systematic uncertainties due to elliptic flow are shown by 
the shaded areas. 
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FIG. 20: (Color online) Fit Gaussian widths of the near-side jet-like correlation, the near-side ridge, and the away-side conical 

emission: (a) as a function of (f> s = \4>t — ^ep| for 1 < pf < 2 GeV/c, and (b) as a function of p^ for the s -integrated 
correlation. The Ar] Gaussian width for the jet-like correlation is also shown. Data are from 20-60% Au+Au collisions. The 
trigger range is 3 < p!j7 < 4 GeV/c. Error bars are statistical only. The near-side jet-like A(f> correlation Gaussion width 
from the minimum bias d+Au data is indicated by the arrow in (a) and by the shaded area in (b); the width of the shaded 
area indicates the statistical uncertainty. 
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FIG. 21: (Color online) Away-side double-peak position relative to A(f> = it from four-Gaussian fits to A(f> correlations at 
A77I > 0.7: (a) as a function of (f> s — \4> t — i^ep\ for two associated particle p^ bins, and (b) as a function of p^ ' for (dis- 
integrated as well as in-plane and out-of-plane correlations. Data are from 20-60% Au+Au collisions. The trigger p^' range is 
3 < py' < 4 GeV/c. Error bars are statistical. The systematic uncertainties due to elliptic flow are indicated by the dashed 
lines. 



obtained by the difference between small and large A77 
correlations by Eq. (TT2"]) ] with a single Gaussian centered 
at A(f> — 0. The fitted Gaussian widths are shown in 
Fig. [201 by the shaded area whose vertical breaths indi- 



cate the statistical uncertainties. The p^-intcgrated cor- 
relation Gaussian width is shown by the arrow in the left 
panel. As seen from the figure, the Acf> widths of the near- 
side jet-like correlations in Au+Au collisions are consis- 
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tent with those from d+Au collisions at the correspond- 
ing pt- In addition, as shown in Fig. [151 the near-side 
jet-like yields are the same between Au+Au and d+Au. 
In fact, the near-side jet- like correlations in Au+Au colli- 
sions of all <j) s bins are consistent with the minimum bias 
d+Au data for all p^ and bins, as shown in Figs. [501 
and This strongly suggests that the near-side jet-like 
correlations in Au+Au collisions are result of in- vacuum 
jet fragmentation, just as in d+Au collisions. 

Figure [21ja) shows the conical emission peak angle as 

a function of <j> s for two associated particle p^ bins. The 
peak angle is approximately constant for <j) s < 45°. For 
4> s > 45° it increases with increasing (f> s , and becomes 
somewhat different for low and high associated particle 
. The larger conical emission peak angle for out-of- 
plane trigger particles may be due to a more significant 
influence from medium flow. For the in-plane orienta- 
tion, the conical emission hadrons from the away side 
are likely aligned with the medium flow, receiving only 
a small deflection to their px- Moreover, the overlap 
collision zone is thinner in the in-plane direction, thus 
the away-side correlated hadrons can escape the collision 
zone more easily. For the out-of-plane orientation, on 
the other hand, the conical emission hadrons move more 
or less perpendicularly to the medium flow direction be- 
cause of the long path-length they have to traverse. They 
receive a large side-kick from the medium flow, broaden- 
ing their final emission angle. 

Figure I2~l7 b) shows the conical emission peak angle as 

a function of p^ for in-plane, out-of-plane, and all trig- 
ger particle orientations. The peak angle is relatively 
independent of the associated particle p^ f° r in-plane 
trigger particles. They may more closely reflect the av- 
erage emission angle of correlated away-side hadrons be- 
cause the medium flow effect is expected to be small, as 
discussed above. 

On the other hand, the emission peak angle for the 
out-of-plane orientation is larger, consistent with a larger 
broadening effect from the medium flow. However, the 
angle position increases with pi"', which is naively not 
expected if those particles are pushed by media with the 
same flow velocity. We note that the medium flow can 
either broaden or shrink the conical emission angle, de- 
pending on the relative orientations of the conical emis- 
sion direction and the direction of the flow. Investigations 
of medium flow effects on the conical emission require re- 
alistic dynamical modeling which is outside the scope of 
this paper. 

It is also worthwhile to note that the peak positions re- 
ported here are from fits to dihadron correlations. They 
are different from those obtained from three-particle cor- 
relations [2l[ where the conical emission angle was found 

to be independent of the associated particle p^\ The 
conical emission angle from the three-particle correlation 
fit is cleaner because the peaks are more cleanly separated 
in the two-dimensional angular space, while the fit to di- 
hadron correlations is more affected by other physics ef- 



fects. One such effect is jet deflection [jl, [7||, which was 
found to be present by three-particle correlations where 
the diagonal peaks are stronger than the off-diagonal con- 
ical emission peaks [21j . 



V. CONCLUSIONS 

Dihadron azimuthal correlations in non-central 20-60% 
Au+Au collisions are reported by the STAR experiment 
as a function of trigger particle azimuthal angle rela- 
tive to the event plane (4> s = \4>t — ^epI) in six equal- 
size slices. The <p a dependence of the dihadron corre- 
lation signal, as well as the trigger and associated par- 
ticle transverse momentum (pt) dependences, are stud- 
ied. Minimum-bias d+Au collision data are presented for 
comparisons. The correlation functions are also obtained 
for small and large pseudo-rapidity separations (|Ary|) in- 
dependently in order to isolate the jet-like and ridge (long 
range A?y correlation) contributions. The resulting jet- 
like and ridge components are studied as a function of 
4> s , trigger particle p^ and associated particle p^ . 

The zyam background subtraction method is described 
in detail. The flow correction is carried out to the order 
of V2V4. The systematic uncertainties in the background 
subtraction are extensively discussed. The effect of possi- 
ble triangular flow fluctuations is not subtracted. The tri- 
angular flow effect, just as the elliptic flow, may be dom- 
inating the untriggered two-particle correlations without 
a high pt trigger particle. However, in our dihadron cor- 
relations with high pt trigger particles, the effect of the 
triangular flow is relatively small. We estimated the tri- 
angular flow effect using both AMPT and UrQMD model 
results as well as the untriggered two-particle correlation 
measurements. The estimated v% maximum-minimum 
modulation is only about 20% compared to our dihadron 
correlation amplitude. 

The dihadron correlations are strongly modified in 
Au+Au collisions with respect to minimum-bias d+Au 
collisions. The modifications strongly depend on the 
trigger particle orientation relative to the event plane 
and evolve with associated particle pft • No significant 
changes are observed between trigger particle p^ ranges 
of 3 < p { T> < 4 GeV/c and 4 < < 6 GeV/c. 

The <j) s and p^ ' dependences of the correlation functions 

are quantitatively similar in the two trigger particle pj) 
ranges. 

The away-side dihadron correlation broadens from in- 
plane to out-of-plane by nearly a factor of 2 in RMS. The 
away-side correlation for <j) s < ir/6 is single-peaked, in- 
dependent of , and not much wider than in d+Au, 
while the amplitude is larger than the d+Au data. For 
4> s > 7r/6, the away-side double-peak structure starts 
to develop and becomes stronger for increasing <f> s and 
increasing p^\ The strongest double-peak structure is 
found at large p^P in the out-of-plane direction. 
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The away-side dihadron correlation amplitude at A<j) = 
tt drops from in-plane to out-of-plane, while that in the 
cone- region remains approximately constant over <j> a . For 
in-plane <f> s , the amplitude ratio in the 7r-region to the 
cone- region increases with p^\ consistent with rf+Au 
and qualitatively consistent with punch-through jets or 
away-side jets not interacting with the medium. How- 
ever, the individual amplitudes in these two regions are 
both higher than in cJ+Au, suggesting other physics 
mechanisms at work. For out-of-plane <f> a , the amplitude 
ratio decreases strongly with p^\ opposite to expecta- 
tions from punch-through jets. 

The trends of the away-side modification underscore 
the importance of the path-length that the away-side par- 
ton transverses in the medium. The away-side medium 
path-length in the reaction-plane direction in 20-60% 
Au+Au collisions is modest and may not be enough to 
generate significant modification to the away-side corre- 
lation. The strongest modification is found for trigger 
particles perpendicular to the reaction plane where the 
away-side medium path- length is the longest. 

The near-side dihadron correlation amplitude de- 
creases with increasing <f) a . The decrease comes entirely 
from the decrease in the ridge. The ridge is extracted 
from correlations at |A»7| > 0.7. Its amplitude is found 
to decrease with increasing (f> s significantly in the 20- 
60% centrality. A ridge with much smaller amplitude 
is found for trigger particles perpendicular to the event 
plane. This feature is present for all associated particle 
p^ , and appears to be independent of pfi . 

The jet-like contribution to the near-side correlation is 
extracted from the difference between small- and large- 
Ar) azimuthal correlations, subject to small experimental 
systematic uncertainties. The jet-like contribution is con- 
stant from in-plane to out-of-plan within our systematic 
uncertainties, and is found to be the same as in d+Au 
collisions. 

The different behaviors of the jet-like component and 
the ridge with respect to (f) s suggest that their production 
mechanisms are different. The jet-like component is in- 
sensitive to the reaction plane and appears to be univer- 
sal, suggesting in-vacuum jet-fragmentation of partons 
whose production is biased towards the surface of the 
collision zone by requiring the high px trigger particles. 
The strong dependence of the ridge on the reaction plane 
suggests its origin to be connected to the medium, not to 
the jet. 

It is found that the near-side jet-like yield and the 
away-side cone yield both have little dependence on <f> s . 
The jet-like spectral shape and conical emission hadron 
spectral shape do not change with <f> s ; however, the cone- 
region spectra are softer than the jet-like spectra. The 
jet-like spectrum for 4 < pj) < 6 GeV/c is somewhat 

harder than that for 3 < p^ < 4 GeV/c, while the 
cone spectra remain the same between the two trigger 
Pj) ranges. These results suggest that the away-side 
conical emission hadrons are related to the away jet part- 



ner, in approximate energy balance with the near-side jet, 
but their production mechanism is different from that of 
jet fragments. The near-side jet-like correlated hadrons 
are likely produced from in-vacuum fragmentation, while 
correlated hadrons in the away-side cone-region likely 
emerge from conical flow, which is related to medium 
properties. 

On the other hand, the near-side ridge and the away- 
side 7r-region appear to trace each other as a function of 
4> s , and also approximately as a function of pft , suggest- 
ing the possibility of a back-to-back ridge. This would 
be consistent with the recent suggestion that the ridge 
may be generated by fluctuations in the initial color flux 
tubes focused by transverse radial flow. Such a picture 
would also explain the decreasing ridge from in-plane to 
out-of-plane because both the color flux tubes and the 
radial flow are the strongest along the in-plane direction. 
However, it remains unclear why the ridge particles are 
much harder than inclusive hadrons. 

The dihadron correlation structure is fit with a four- 
Gaussian model representing a back-to-back ridge and an 
away-side conical emission. The fitted conical emission 
angle increases from in-plane to out-of-plane. The larger 
conical emission angle out-of-plane may be due to the ef- 
fect of medium flow. For in-plane trigger particles, the 
away-side conical emission angle is approximately con- 
stant over the associated particle p^ . For out-of-plane 
trigger particles, the away-side conical emission angle in- 
creases with pj^ . The larger emission angle at large Pj^ 
seems counterintuitive as a result of broadening by the 
medium flow. More studies are needed to further our 
understanding of the away-side correlation structure. 

To summarize our main findings, highpr triggered jets 
are biased towards surface emission, and the jet fragmen- 
tation is hardly modified by the medium. Away side part- 
ner jets interact maximally with the medium in the di- 
rection perpendicular to the reaction plane. These inter- 
actions result in a conical emission of correlated hadrons 
and the strong double-peak correlation structure on the 
away side. The near-side jet-like component and the 
conical emission hadrons appear to have a fixed rela- 
tive energy balance. The near-side jet-like component 
is accompanied by the ridge in the reaction-plane direc- 
tion. The ridge magnitude drops rapidly with increas- 
ing <f> s and largely disappears out-of-plane in mid-central 
20-60% Au+Au collisions. The away-side 7r-region am- 
plitude appears to trace the near-side ridge amplitude, 
suggesting a back-to-back ridge. The large hadron yield 
in the 7r-rcgion with in-plane trigger particles results in 
a single correlation peak on the away side, even though 
the cone hadron yield is as appreciable as for other trig- 
ger particle orientations. The most natural explanation 
for our results seems to be the combination of a near- 
side in-vacuum jet-fragmentation, a back-to-back ridge 
in the reaction plane, and an away-side conical emission 
induced by jet-medium interactions. 
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Appendix A: Dihadron Correlation Functions 

This appendix presents dihadron correlation functions. 
Figures [HI [231 HH and 1251 show the raw correlation func- 
tions. Figures ESI EH Ell and [29] show the flow and 
ZYAM background-subtracted correlation functions. Fig- 
ures [30] and [31] show the near-side jet-like correlation 
functions. The data for the correlation functions and 
all other figures in the paper are published on-line at 
|http://www.star.bnl.gov/central/publications/| 
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FIG. 22: (Color online) Raw dihadron A<j> correlations with trigger particles in six slices of azimuthal angle relative to the event 
plane, 4> s = \ 4>t — ipEP j . The data are from minimum-bias 20-60% Au+Au collisions. The trigger pr range is 3 < p^} < 4 GeV/c. 

Five associated particle p£ bins are shown. Both the trigger and associated particles are restricted within |r/| < 1. Statistical 
errors are smaller than the symbol size. The curves are flow modulated zyam background calculated with four-particle «2{4} 
and away-side two-particle «2{2,AS} (dashed curves), and the average V2 from the two methods (solid curve). 
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FIG. 23: (Color online) Same as in Fig. [22] but for trigger particle 4 < < 6 GeV/c and six bins in associated particle 
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FIG. 24: (Color online) Same as in Fig.[22]but for |Ar?| > 0.7. 
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FIG. 26: (Color online) Background-subtracted dihadron correlations with trigger particle in six slices of azimuthal angle 
relative to the event plane, <j> s — \<pt — V"ep|- The trigger pr range is 3 < < 4 GeV/c. Five associated particle ' bins 
are shown. Both the trigger and associated particles are restricted to be within \r)\ < 1. The figure corresponds to the raw 
correlations in Fig. 1221 The data points are from minimum-bias 20-60% Au+Au collisions. Systematic uncertainties are shown 
in the thin histograms embracing the shaded area due to flow subtraction and in the horizontal shaded band around zero due to 
ZYAM background normalization. Statistical errors are mostly smaller than symbol size. For comparison, the inclusive dihadron 
correlations from d+Au collisions are superimposed as the thick (green) histograms (only statistical errors are depicted). 




FIG. 27: (Color online) Same as in Fig. [26]but for trigger particle 4 < < 6 GeV/c and six bins in associated particle 
The figure corresponds to raw correlations in Fig. 1231 
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FIG. 28: (Color online) Same as in Fig, 1261 but for \Ar]\ > 0.7. The figure corresponds to raw correlations in Fig. 1241 




FIG. 29: (Color online) Same as in Fig, 1271 but for \Ar]\ > 0.7. The figure corresponds to raw correlations in Fig. 1251 
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FIG. 30: (Color online) Jet-like dihadron correlations with trigger particle in six slices of azimuth relative to the event plane, 
4>s = \4>t — V'epI- The jet-like dihadron correlations are obtained from the difference between \Arj\ < 0.7 and (acceptance 
weighted) \Ar]\ > 0.7 correlations. The trigger pt range is 3 < p^) < 4 GeV/c. Five associated particle p£ bins are shown. 
Both the trigger and associated particles are restricted to be within < 1. The data points are from minimum-bias 20-60% 
Au+Au collisions. Superimposed for comparison in the thick histograms are the inclusive jet-like dihadron correlation from 
ci+Au collisions. Errors bars are statistical; Systematic uncertainties are small. 
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FIG. 31: (Color online) Same as in Fig. [30] but for trigger particle 4 < < 6 GeV/c and six bins in associated particle 
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